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ABSTRACT 

Context. During the formation of a star, material is ejected along powerful jets that impact the ambient material. This outflow phe- 
nomenon plays an important role in the regulation of star formation. Understanding the associated shocks and their energetic effects 
is therefore essential to the study of star formation. 

Aims. We present comparisons of shock models with observations of FL and SiO emission in the bipolar outflow BHR71, and predict 
water emission, under the basic assumption that the emission regions of the considered species coincide, at the resolution of currently 
available observations. 

Methods. New SiO observations from APEX are presented, and combined with Spitzer and ground-based observations of H 2 to con- 
strain shock models. The shock regions associated with targeted positions in the molecular outflow are studied by means of a ID code 
that generates models of the propagation of stationary shock waves, and approximations to non-stationary ones. 
Results. The SiO emission in the inner part of the outflow is concentrated near the apex of the corresponding bow-shock that is also 
seen in the pure rotational transitions of H 2 . Simultaneous modelling is possible for H 2 and SiO and leads to constraints on the silicon 
pre-shock distribution on the grain mantles and/or cores. The best-fitting models are found to be of the non-stationary type, but the 
degeneracy of the solutions is still large. Pre-shock densities of 10 4 and 10 5 cm -3 are investigated, and the associated best-model 
candidates have rather low velocity (respectively, 20-30 and 10-15 km s ) and are not older than 1000 years. We provide emission 
predictions for water, focusing on the brightest transitions, to be observed with the PACS and HIFI instruments of the Herschel 
Telescope. 

Key words. Stars: formation - ISM: jets and outflows - ISM: individual objects: BHR71 - Submillimeter: ISM - Infrared: ISM 
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1. Introduction into the gas phase, such as silicon-bearing species (Schilke et al. 

\fe ! ll997HGusdorf et alj|2008al hereafter G08a). The time-scales in- 

m ' The propagation of shock waves is an ubiquitous phenomenon volved in the heatin S and j? some , of th f s } ock chemica l pro- 
OO in the interstellar medium, where the sound speed is low due c f se * £fJ hort (10 2 to 10 4 years; | Gusdorf et al. ||2008b, here- 
O to the temperature conditions. Observations over the past few after G08b )' s ° the shocked region rapidly acquires a chemical 
t-H ; decades ( Snell et al. 1980) have shown that, at the early stages composition distinct from that of the quiescent medium. As the 
, f the formation of solar-type stars, the process of mass ac- 8^ cools down radiatively, reactions with high energy barriers 
cretion is almost always associated with mass ejection, in the no lon g er °P erate ! and * e chemistry is dominated again by low- 
form of collimated jets, extending from a few astronomical units temperature reactions. The resulting molecular emission can be 
up to parsecs from the exciting source (Dougados et al. 2000; used as a dla g n ° stl c ^1 to study the physical and chemical char- 
£ ■ |BallvetalJfT996h . The supersonic impact between the jet and actenstics of the shocked region. 

the parent cloud generates a shock front, which propagates in Recently, G08b placed constraints on shock parameters us- 
the collapsing interstellar gas, and also an inverse shock that ing H2, SiO and CO observations in the Bl knot of the LI 157 
propagates along the jet itself. Large cavities, called bipolar outflow. At the same time, the HIFI and PACS receivers on 
outflows, appear, which have been extensi vely studied thro ugh Hersc hel have observed the entire LI 157 outflow dNisini et alj 
the molecular emission that they generate (Bac hille3ll996h . At 1201 0|) and provided velocity-resolved da ta of the LI 157-B1 knot 
the apices of these cavities, the shock wave heats, compresses dLefloch et al.ll20iot ICodella et alJl2010h . Their studies indicate 
and accelerates the ambient interstellar gas. As the temperature that the H2, SiO and H2O emission regions are spatially coin- 
rises to at least a few thousands degrees, the energy barriers cident on large scales. They also shed light on the necessity to 
of numerous reactions involving neutral or ionized molecules provide reliable models of water emission in order to better un- 
can be overcome an d the chemistry of c ertain species can be derstand its product ion in shocks and its role in t he coo ling of the 
significantly altered (Bachiller et al. 2001). Similarly, processes outflow. Recently. iFlower & Pineau des Foretsl (12010) (hereafter 
specific t o the propagation of shock waves affect the interstel- FP10) presented an extensive study of water emission for a small 
lar gr ains (Fl ower & Pineau des Forets 2003*1; iGuillet et ai]|2007l, grid of stationary shock models, and compared it to Spitzer ob- 
2009), leading to the injection of molecular and atomic species servations of water emission in NGC2071, with a certain level of 
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agreement. Our aim is to push the confrontation between models 
and observations one step further, by using the models, which 
have already demonstrated their ability to reproduce some ob- 
servations, to predict water line intensities. Using an approach 
to comparing models with observations that is similar to that of 
G08b, as applied to H2 and SiO, we aim to constrain the physi- 
cal and chemical parameters of the targeted region. In a second 
step, we incorporate these preliminary constraints and the lat- 
est compilation of molecular data for H2O in an LVG radiative 
transfer code to provide predictions of water emission, under the 
assumption that the emission regions of the various species co- 
incide. 

The choice of target is of crucial importance, as our candi- 
date must have been the object of prior studies, and must be a 
good potential target for contemporary and future instruments 
such as Herschel and ALMA. Located at a dista nce of about 175 
pc, the Bok globule BHR71 dBourke et all 19951) is a well-known 
example of isolated star formation observable in the south- 
ern hemisphere, and hence is a reasonable candidate for future 
ALMA observations. Its observation by the Herschel telescope 
is guaranteed through the key-program WISH, and also DIGIT. 
The associated outflow region contains signatures of two distinct 
molecular jets dBourke et all Il997t iMvers & Mardonesl 1 1998; 
iParise et alJ2006|) centered on two different protostellar sources, 
IRS1 and IRS2 (|Bourkej|20011) separated by -3400 AU. Bright 
HH objects are associated w ith the blue-shifted lobes o f each 
outflow, HH320 and HH321 (ICorporon & Reipurthll 19971) . They 
have been imaged in the S n transition at 67 1 1 A, indicating that 
at least part of the outflows are dissociati ve. The dynamical age 
of HH321 is estimated to be 400 years (ICorporon & Reipurthl 
119971) . It has not yet been possible to determine the dynamical 
age for HH320. NIR JHK-ban d spectra cove r ing bo th HH320 
and HH321 were obtained by Giannini et al. (2004). Through 
detailed shock modelling they derived a pre-shock density of 
10 4 crrT 3 and shock velocity of 41 km/s for HH320A; this was 
found by fitting a non-steady-state shock model to the observed 
H2 lines. The age of the non-steady-state shock is 475 years, 
which is more or less in agreement with the dynamical age of 
the HH321 flow. The pre-shock density is lo wer than is predicte d 
on the basis of CO observations, 10 cm " 3 (IParise et al.l l2006). 
However, the latter is the density of the molecular outflow, which 
is compressed relative to the ambient pre-shock cloud. These in- 
ner parts of the bipolar outflow, covering the HH objects, have 
also been targeted by the Spitzer Space Telescope, and the cor- 
responding results in terms o f rotational H? excit ation diagrams 
have been made available by iNeufeld et alj (120091) . Another ar- 
ticle is dedicated to the p arameters map analysis base d on these 
data dGianninietal .1120111) . Finally. iGarav et al.l (1 19981) have pre- 
sented maps of SiO (2-1) and (3-2) of the whole bipolar outflow. 

In Section|2] we presented a c ombined analy s is of b oth rovi- 
brational transitions of H2 from iGiannini et al.l (120041) an d ro- 
tational lines, as compiled by Spitzer dNeufeld et al.l l2009). We 
then briefly present the code that is used to simulate the propaga- 
tion of shock waves through the interstellar medium, and intro- 
duce the grid of models that we computed. Finally we describe 
the constraints that we obtain in several knots of the BHR7 1 out- 
flow by comparing more or less complete H2 observations with 
models. In Section [3] we present the results of comparisons be- 
tween observations of SiO and the models. We describe the ob- 
serving setups, the actual observations, and the results in terms 
of simultaneous fits of the H2 and SiO molecules. Section [4] 
contains our predictions of water emission line intensities. We 
briefly summarize our methodology and then discuss the appro- 
priateness of the steady-state and statistical equilibrium assump- 



Table 1. The considered positions in the BHR71 outflow. 
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tions, before presenting a comparison of results obtained using 
each of these approximations. Section[3]contains our concluding 
remarks. 

2. H 2 observations and modelling 

2.1. Observations 

There are two distinct sets of H2 observations. The details of the 
ground-b ased observations of th e rovibrational transitions are re- 
ported in lGiannini et all (I2004I) . The Spitzer observations of the 
pure rotational transitions, described in Neufeld et al. (2009), are 
the subject of an separate article (Giann ini et alJ[201lb . 

The comparison between observations of H2 emission and 
model results is usually based on excitation diagrams that are 
derived for selected emission regions. Such excitation diagrams 
display ln(N V j/gj) as a function of E V j/k^, where N V j (cm -2 ) is 
the column density of the rovibrational level (v, J), E V j/kn is its 
excitation energy (in K), and gj = (2j + 1)(2I + 1) its statistical 
weight (with 1=1 and / = in the respective cases of ortho- 
and para-H2). If the gas is thermalized at a single temperature, 
all the points in the diagram fall on a straight line. 

In their original s tudy of conditions in the BHR71 outflow, 
Gia nnini et al. I d200l extracted the excitation diagrams corre- 
sponding to nine positions (knots 1 to 9) and four HH objects 
(HH320A/B and HH321A/B), in order to compare their obser- 
vations with shock models. In the present study, we attempt to 
make use of most of the available data, both in the pure rota- 
tional and in the rovibrational transitions of molecular hydro- 
gen. Unfortunately, combining the two sets of observations is 
only possible where both sets of data are available and when 
a cross calibration can be performed to ensure their consistency. 
We consequently make use of the following excitation diagrams : 

- pure rot ational transitions for the 'H2 knots' 5, 6, 8 as de- 
noted by IGiannini et all d2004l) ; 

- rotational and rovibrational transitions for the two Herbig- 
Haro obje cts that are cal l ed HH 320AB and HH321AB, as 
defined by IGiannini etal.1 d2011l) ; 

- pure rotational transitions on the 'SiO knot' (see 
Subsection l3~2i >. 

The corresponding list of coordinates is provided in TableQ] 
All these points are presented in FigureQ] which shows the inter- 
nal parts of the BHR71 outflow, as mapped by the Spitzer Space 
Teles cope in the 0-0 S(5) spectral line of H2 (Neufeld et al. 
2009). In each case, the excitation diagrams were extracted for 
a circular emission region of 5" radius, centered on these po- 
sitions, shown for the case of the so-called 'SiO peak' (see 
Subsection 13. 2\ in Figure Q] (smallest circle). In the case of the 
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Fig. 1. The internal parts of the BHR7 1 outflow as observ ed by 
Spitzer in the H 2 0-0 S(5) transition (Ne ufeld et al.ll2009l) : the 
colour scale unit is erg cirT 2 s~' sr _1 . The contours correspond 
to the levels of 3, 9, 15, 21, 27, and 33 cr. The pink hexagons give 
the positions of the IRS1 and IRS2 sources as determined by 
iBourke et al.l (2001). The positions of the Herbig-Haro objects 
HH320AB and HH321AB are plotted in black as well as those 
of the SiO peak (referring to the point of maximum SiO (5-4) 
emission, see Section|3]l, and of the SiO knot, which was chosen 
to perform our simultaneous fit of H2 and SiO emission, and 
also to make predictions of water emission. The other emission 
knots that were analysed are plotted in white (k nots 5, 6, and 
8, following the labelling of iGiannini et al.l [20041) . Around the 
SiO-peak, four radii are displayed, corresponding to 5, 9, 12, 
and 14.4". The first represents the size of our region of analysis. 
The latter three correspond to the beam of the APEX telescope 
for the SiO (8-7), (6-5) and (5-4) transitions, respectively; see 
text. 



Herbig-Haro objects, the consistency between the rotational and 
rovibrational data sets was ensured by means of the following 
method, as explained by IGiannini et al.l d201 ll) . First, the image 
of H2 1-0 S(l) line at 2.122 pm (retrieved from the ESO archive 
at http://archive.eso.org/eso/esojrchive_main.html) is degraded 
to the Spitzer spatial resolution of 2"/px. Then, this image is cal- 
ibrated by using 2MASS fluxes of several bright stars in the field, 
and the photometry is computed in the same areas of 5" radius 
as the Spitzer images. To find the intercalibration factor between 
imaging and spectroscopy, the measured 2.12 pm flux was corn- 



Table 2. Input parameters of our grid of stationary C- and non- 
stationary CJ-type shock models. 
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pared with the flux measured in the 1" slit where spectroscopic 
observations have been performed. Finally, the same factor was 
applied to the fluxes of the other NIR lines, adopting averages of 
the line ratios (relative to the 2.12 pm transition) in the A and B 
substructures of HH320 and HH321. 

The data were corrected for extinction, adopting the upper 
limit of the visual extinction from lGiannini et a l]( 2004). A y _ 2, 
and u sing the interstellar extinction law of Rieke & Lebofsky 
(1985) . We note that the effect of the interstellar extinction on the 
observational points is minor. The excitation diagrams resulting 
from the observations are plotted in Figure|2] The rotational and 
rovibrational temperatures derived from these observations are 
also indicated on the plots in the same colour code, along with 
the corresponding line. Again, the reddening correction has only 
minor significance on the temperature values. As anticipated, 
the rovibrational transitions are found to trace warmer gas than 
the pure rotational transitio ns (see, for example. lLe Bourlot et al.l 
12002- IGiannini et alJl2006l and G08b). 

2.2. Grid of models 

The shock code that we use to simulate the propagation of shock 
waves in the interstellar medium is described in G08b. It gen- 
erates one-dimensional stationary models of both continuous 
or discontinuous kind (respectively, C- or J- type shocks), and 
also approximations to time-dependent solutions, in the form 
of a shock wave p ossessing the characteristics of both C- a nd 
J-type models (see IChieze etai1ll998r. iLesaffre et al.l l2004alTbh . 
The shock code solves the differential equations that determine 
the evolution of a set of dynamical and chemical variables: the 
neutral and ionized fluid temperatures, velocities and densities; 
the fractional abundances of over 125 species, linked by more 
than 1000 chemical reactions. The populations of the first 150 
levels of molecular hydrogen ar e calculated with i n the shock 
code, based on the treatment of iLe Bourlot et all (120021) . with 
the only difference being that we adopted the most r ecent H- 
H2 collisional rate coefficients (Wrathmall et al. 2007). Finally, 
grain-related processes are also in cluded in the code, following 
iFlower & Pineau des Foritsl (120031) . 

The parameters that we aim to constrain are the inputs of the 
models, namely: 

- ambient-medium-related quantities: pre-shock den- 
sity «h and magnetic-field parameter b, defined by 
B(pG) = b x V«h (cm- 3 ); 

- shock-related quantities: shock velocity v s , and age in the 
case of a non-stationary solution; 

- chemistry-related quantities, such as the initial value of the 
ortho-para ratio for molecular hydrogen, or the initial distri- 
bution of silicon-bearing material in the grains. 

As a first step, we made use of a grid of stationary C- 
type shock models with the parameters specified in Table [2] We 
should mention that the parameter coverage is not complete, in 
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Fig. 2. The H2 excitation diagrams observed for each position indicated in Figure Q] in blac k (pure rotational li nes), and grey 
(rovibrational lines, when available). The data were corrected for extinction, using Ay = 2 from lGiannini et aLl (I2004I) . The inferred 
rotational temperatures are indicated in black, and the rovibrational temperatures, in grey. The best-fitting model results are also 
shown in each case, in blue and/or red symbols. The shock parameters are given in each panel in the corresponding colour. 



the sense that not all velocities are present in our grid for all val- 
ues of the magnetic-field parameter, b. In fact, the velocity of 
C-type shocks must remain below a critical value that depends 
main ly on the pre-shock density and magnetic-field parame - 
ter dLeBourlot et al.ll2002t iFlower & Pineau des Foritsl I2003I) . 
which explains the decrease of the maximum shock velocity with 
the pre-shock density in our grid. Following the method pre- 
sented in G08b (Section 4.1), this set of C-type shock models 
enabled us to restrict the range of the search in the parameter 
space for the CJ-type shock models. We then computed a grid of 
non-stationary shock models around a first estimate of the shock 
age, making sure that the range of ages was sufficient to include 
any model likely to fit the H2 observational data. The grid of 
CJ-type models covers the combination of parameters given in 
Table |2 for ages ranging from a few tens to a few thousands of 
years. 

For each of the above models, we obtained an excitation di- 
agram that can be compared with observations. 

2.3. Comparisons and initial constraints 

We describe here the results of the comparisons between obser- 
vations and our models. A few preliminary remarks about our 
results are applicable to all the knots and regions analysed. First, 
the filling factor was considered to be equal to 1 for all transi- 
tions; this yields H2 knot sizes consistent with the thickness of 



the H2 emitting layer of our models. Secondly, under this as- 
sumption, we find that, in all cases, only non-stationary shock 
models reproduce satisfactorily both the pure rotational and rovi- 
brational (when available) parts of the H2 excitation diagram, 
a result that has already been established in similar studies of 
other objects (e.g. HH54. lGiannini et al.l2006l orL1157-Bl, see 
G08b). Finally, the departure of the observational points from a 
straight line in the purely rotational parts of the excitation di- 
agram ('saw-tooth' pattern, which is weak in our cases) is in- 
dicative of deviations of the ortho-to-para ratio from the high- 
temperature LTE value of 3.0. Such a pattern can be reproduced 
in the models by slightly modifying the value of the initial ortho- 
to-para ratio; but fine-tuning this parameter is not a major con- 
sideration in the context of the present study. 

We give prior focus to the Herbig-Haro objects, for which 
we can benefit from H2 excitation diagrams combining pure ro- 
tational and rovibrational transitions, and to the so-called 'SiO 
knot', for which future observations of SiO rotational transitions 
are expected to provide further constraints on our shock mod- 
els. Tables lATI and IA~2l present a summary of the restrictions 
that comparisons between observational and computed excita- 
tion diagrams place on the shock model parameters. In these 
tables, the values of the pre-shock density, magnetic-field pa- 
rameter and shock velocity are indicated for our best models, 
as well as the age range for which the excitation diagrams are 
in satisfactory lagreement with the observations. In the case of 
the Herbig-Haro objects HH320A/B and HH321A/B, the best- 
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fitting models are non-stationary, and the pre-shock density is 
constrained to 10 4 cirT 3 , which is in agreement with the average 
density of the who le globule (9xl0 3 cm 3 ), as determined by 
Bour keetalJdl997h . Shock velocities in the range 15-35 km s 
yield good fits to the observations, depending on the associated 
magnetic-field parameter and age. We note that the shock ages 
that we obtain are of the same order of magnitude as the observa- 
tional dynamical ages. The excitation conditions are remarkably 
similar in both sources, in spite of their location on a different 
lobe of the outflow, and the probable mixing of the two outflows 
from IRS1 and IRS2. In the case of HH320A/B, o ur findings 
are in agreement with those of Gi annini et al.l d2004) in terms of 
pre-shock density, magnetic field, and shock type and age, but 
our value of the shock velocity is smaller than theirs. Several 
factors could account for this difference, which might not be a 
real discrepancy, given the degeneracy of our fitting models. The 
first factor is that their value lies slightly outside the range of 
our parameter search. In addition, we studied HH320A/B, rather 
than HH320A. Furthermore, we added the pure rotational tran- 
sitions of H2 to the analysis. Finally, our shock model incorpo- 
rat es an update of th e H-H 2 collisional rate coefficients, follow- 
ing |WrathmaiLetaD (120071) . 

The parameters of the SiO knot are not as well constrained, 
as only the pure rotational observations are available to build the 
H2 excitation diagram. Similarly to the HH objects, the best fit- 
ting models are of non-stationary shocks with velocities in the 
range of 15-35 km s _1 for a pre-shock density of 10 4 cirT 3 (with 
ages dependent on the value of the magnetic-field parameter). 
Additionally, a few models with a pre-shock density of 10 5 cm -3 
are found to fit the observations equally satisfactorily, with even 
more moderate shock velocities (10-20 km s ) and lower shock 
ages (which again depend on the value of the magnetic-field pa- 
rameter). Examples of the comparisons between observations 
and models are shown in the three upper panels of Figure |2 
for the cases of the SiO knot, HH320A/B and HH321A/B, re- 
spectively. In each panel, the parameters corresponding to the 
selected shock model are also indicated. 

For the three other points studied( 'knot 5', 'knot 6', and 
'knot 8'), once again the excitation diagrams do not constrain the 
parameters as tightly as in the HH regions, since only the pure 
rotational transitions of H2 can be used. The number of models 
that fit these data is large, but we show some examples of the 
best fits that we could derive in the lower panels of Figure|2] We 
note that the modelling work conducted for these positions also 
suggests similar excitation conditions, regardless of their associ- 
ation with one lobe of the outflow or the other. 

3. SiO observations and modelling 

3.1. Observations 

The observations were carried out mainly in two observational 
runs, in August and Oct ober 2009, with th e Atacama Pathfinder 
Experiment (APEXfl iGusten et all 120061) . We used the facil- 
ity's own APEX-1 and APEX-2 receive rs (see respectively 
IVassilev et al.l 120081 iRisacher et al.l 120061) to observe respec- 
tively the SiO (5-4) and (6-5) transitions, and the SiO (8-7) tran- 
sition, in combination with the MPIfR Fast Fouri er Transform 
Spectrometer backend (FFTS, iKlein et all 120061) . Focus was 
checked at the beginning of each observing session, and after 

1 This publication is partly based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 



Table 3. Observed lines, frequencies, and corresponding tele- 
scope parameters: beam and forward efficiencies, beam sizes, 
typical system temperatures, and spectral resolution. 



line 


SiO (5-4) 


SiO (6-5) 


SiO (8-7) 


v (GHz) 


217.105 


260.518 


347.331 


B e fl 


0.75 


0.75 


0.73 


-fcff 


0.95 


0.95 


0.95 


beam (") 


28.7 


24.0 


18.0 


T sys (K) 


155-175 


270-300 


275-335 


Av (km s~') 


2.02 


1.69 


1.26 



sunrise on Mars. Line pointing was locally checked on 07454- 
71 12 or IRAS 15194-51 15. The pointing accuracy was found to 
be of the order of 5" rms. Table [3] contains the main character- 
istics of the telescope for the observed transitions: beam sizes, 
forward and beam efficiencies, typical system temperatures, and 
the spectral resolution that was used for each transition. The 
observations were per formed in position- switching mode using 
the APECS software dMuders et alj|2006l) . with a reference po- 
sition at offset (-600",600") from the central position indicated 
in TableQ] The data were reduced with the CLASS software (see 
|http://www.iram.fr/I RAMFR/ GILDAS) . 

3.2. Results 

The observations resulted in one SiO (5^4) map of the inner, 
upper (red) lobe, as well as in SiO (6-5) and (8-7) spectra, ob- 
tained by pointing the telescope at designated positions inside 
this lobe. The absolute coordinates of the (0",0") position are 
given in Table Q] 

The map of the SiO (5-4) line covers a roughly rectangular 
region, extending from the lower left offset (35", 25") to the up- 
per right (-60", 170"), and was performed in raster mode, spend- 
ing 30 seconds on each point. The beam sampling was 10", 
which is less than the half beam size of 14.35", meaning the 
map is fully Nyquist sampled. The result is plotted in Figure [3] 
the SiO emission, integrated between -10 and 60 km s" 1 , is dis- 
played in dark blue contours overlaid on the H2 0-0 S(5) map 
already shown in Figure Q] The half-maximum signal contour is 
shown in red. The r.m.s. of the SiO (5-4) map was found to be 
0.26 K km s _I at the resolution 2.02 km s _1 , and the contours 
correspond to levels of 6 to 30 cr, in steps of 3 <x. The contours 
show that the SiO (5^1) emission peaks outside of the Spitzer 
spectral line maps, at a position that will be referred to as the 
'SiO peak', with offset (-15", 87"). In order to perform a simul- 
taneous study of the H2 and the SiO emission, the point to be 
studied must be located inside the Spitzer map: we chose the 
point of maximum emission of H2, at an offset of (0", 77"), to 
which we refer as the 'SiO knot'. We consider the shift between 
the APEX SiO (5-4) and Spitzer H 2 0-0 S(5) emission peaks to 
be real, as it is far too large to be attributable to pointing errors. 
Indeed, the distance between the south-eastern and the north- 
western tips of the half-maximum SiO emission (roughly of re- 
spective coordinates (0",70") and (-30", 110"), see Figure© is 
about 50", which amounts to 40" once deconvolved from the 
SiO (5^4) telescope beam size. Consequently, we consider that, 
in this direction, which coincides with the line between our SiO 
knot and peak (see Figure [3), the emitting region is extended 
enough to allow us to consider that the signal observed at the 
position of the SiO knot is free of any residual signal from the 
SiO peak. Moreover, from the extension of the Spitzer maps, 
it is impossible to exclude the presence of another, potentially 
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Fig. 3. Upper internal part of the BHR71 outflow. The back- 
ground map is the same as in Figure Q] The SiO (5-4) contours, 
which are overlaid in a dark blue colour, correspond to levels 
of 6 to 30 cr, in steps of 3 cr, for an integration between -10 
and 60 km s _1 . The half-maximum, 15 cr contour is displayed 
in red. The associated beam size is indicated in the upper left 
corner. Our SiO (5-4) coverage is actually wider than the field 
displayed on this Figure. 



stronger peak of H2 emission, coinciding with the SiO peak. In 
order to check, we overlaid our map on the four (way more ex- 
tended) channel maps of IRAC, especially the fourth, at 8 mi- 
crons, whose broadband intensity is dominated by the 0-0 S(4 ) 
and S(5) transitions of H2 in this region (lNeufeldetalj r2009). 
We found that these overlays fail to confirm the existence of an 
H2 peak coinciding with that of the SiO (5^) emission. 

The spectra that we obtained of the SiO (5^4) transition at 
both the SiO peak and knot are presented in Figure|4] along with 
the spectra of the other two SiO transitions. In the case of the 
SiO (5-4) transition, the original beam size was 28.7". In order 
to derive the spectra at the two selected positions as exactly as 
possible, we resampled the data with a spatial resolution of 5" 
and then averaged the four spectra closest to the two positions 
to obtain the results that are presented on Figure [4] In the case 
of the SiO (6-5) and (8-7) transitions, in addition to this spatial 
resolution problem, we had also to convolve the data to the SiO 
(5-4) beam size in order to render the whole set of transitions 
consistent. After this preliminary convolution, we obtained the 
spectra displayed in Figure [4] Our mini-maps around the two 
positions were extended enough to cover the larger beam of the 
SiO (5^1) observations. 

The shape of the line profiles is typical of bipolar out- 
flo ws (similar, for in stance, to those acquired around LI 157: 
see Nisini et al .1 120071) . peaking near the cloud velocity and with 
the main signature being the line wings extending to about 
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Fig. 4. The spectra obtained in the three observed SiO transitions 
at the positions of the SiO peak and knot, of respective offsets (- 
15", 87") and (0", 77"). Only the SiO knot is considered in our 
modelling. For each spectrum, we adopted the SiO (5^) spatial 
resolution. The ambient cloud velocity, -4.5 km s _I is indicated 
by the vertical dashed lines. 



50 km s _1 , most visible in the SiO(5-4) line profile of the peak. 
Such wings seem to be present, though narrower, at higher fre- 
quencies. However, it is then more difficult to estimate their 
width, as the signal-to-noise ratio decreases as the observed fre- 
quency increases. The observed SiO probably arises from a com- 
plex mixture of shocks with different velocity projections on the 
line of sight. The velocity range of our shock models does not 
exceed the terminal velocity of 50 km s _I . On the SiO (5-4) and 
(6-5) spectra associated to the SiO peak, a velocity bullet seems 
to arise between 20 and 50 km s _1 . It could arise from a process 
that alters the SiO emission within the region of the shock con- 
sidered, such as an abundance effect, or from the presence of an- 
other shock, with a different velocity, within the telescope beam. 
These possibilities are discussed in Subsection l3.3l but cannot be 
thoroughly investigated, given the large beam sizes associated 
with these transitions; they call for a more detailed study, with 
higher spatial resolution. 

The last step was to derive integrated intensities from these 
spectra. The results for the SiO knot, which is considered in our 
subsequent analysis, are shown in Table [4) integrated intensities 
and associated r.m.s. are given. The velocity interval over which 
the integration was made was [-10, 60] km s for all three tran- 
sitions. Again, the values inferred are co nsistent with resul ts for 
other outflows, such as L1448 or LI 157 dNisini et al.ll2007l) . for 
which the central source is also a young Class proto-star with a 
similar bolometric luminosity (respectively, aroun d 13.5, 11, and 
1.5-3 .5 L for BHR71, LI 157, and L1448; see iNeufeld etail 
2009). In particular, the line profiles we observe for SiO (5^1), 
(6-5), and SiO(8-7) are remarkably similar to those of LI 157- 
Bl, peaking near the cloud velocity. The integrated intensities 
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Table 4. Integrated intensities and r.m.s. for the three observed 
SiO transitions at the position of the SiO knot, of relative coor- 
dinates (0", 77"). 



line 


SiO (5-4) 


SiO (6-5) 


SiO (8-7) 


J 7mbAd (K km s _1 ) 


1.1 


4.1 


1.8 


"7 Jim A, (Kkms- 1 ) 


0.4 


0.4 


0.2 



we inferred are also very clos e to those compile d for each knot 
of the LI 157 bipolar outflow (Nisi ni et al1l2007h . 



3.3. Modelling results 

We derive emission-related quantities by means of a radiative 
transfer code used in combination with the shock models pre- 
sented above. The radiative transfer code is based on the LVG 
approximation, and has been described and used in G08a and 
G08b. At each point of the shock model, it uses the computed pa- 
rameters of relevance, such as the temperature, density, velocity 
gradient, and collision-partner density, to solve the equations of 
statistical equilibrium and calculate emission-related quantities, 
such as the level populations, optical depths, and local emissivi- 
ties. The set of equations is solved by means of a simple lambda- 
iteration, that stops when the relative difference between each 
of the level populations computed in two successive iterations 
drops below a convergence value of 10 4 . The set of collisional 
rate coefficients for SiO with H2 as a collision partner is simply 
scaled by the square root of the re duced mass ratio from the SiO- 
He rate coefficients calculated by Dayou & Balanca (2006). The 
LVG code also provides the integrated intensity over the whole 
shock model for all transitions considered. 

The silicon chemistry is extensively discussed in G08a and 
G08b. Regarding the initial distribution of silicon-bearing mate- 
rial, G08a made the assumption that Si is initially present exclu- 
sively in the grain cores, in the form of silicates, with a fractional 
abundance of 3.35 x 10 relative to total H. G08b introduced an 
alternative scenario, in which the Si is partly in the grain man- 
tles, in the form of SiO. Such an assumption was found to be 
necessary to fit both the SiO and H2 line intensities by means of 
a single shock wave model. Furthermore, a simultaneous fit can 
be achieved only by means of a non-stationary, CJ-type shock. 
When computing CJ-type models in the present study, we adopt 
one of three values of the fraction of elemental silicon initially 
in the form of SiO in the grain mantles, namely 0.00, 0.01, or 
0.10. 

In the first case, when no SiO is initially present in the grain 
mantles, the gas-phase silicon is produced by the erosion of the 
grain cores only, which occurs, in the C- and CJ-shock mod- 
els, where the predominantly charged grains collide with the 
neutral particles at the ion-neutral drift speed (ambipolar diffu- 
sion). Three reaction paths then lead to the formation of SiO in 
the gas phase, namely Si + (OH,H)SiO + (H 2 ,H)HSiO + (e,H)SiO, 
Si(OH,H)SiO, and Si(0 2 ,0)SiO; the last two reactions pre- 
dominate. Following G08a, we adopted the rate coefficient 
that was exper i menta lly measured for the Si+02 reaction by 
iLe Picard et all d2001l) for both Si+0 2 and Si+OH. In the sec- 
ond case, in which SiO is initially present in the grain mantles, 
their rapid erosion by sputtering in the shock wave generates SiO 
directly in the gas phase (see also the discussion in G08b). 

O2 is a key reaction partner in the process of formation of 
SiO in the gas phase, and so its abundance must be modelled to 
the greatest possible accuracy. Furthermore, the question of the 
pre-shock abundance of O2 in the gas phase must be considered 
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Fig. 5. Examples of fits of the integrated intensities of our three 
transitions of SiO. On the y-axis is the integrated intensity 
(K km s~'), and on the x-axis is the rotational quantum num- 
ber of the upper level of the transition. The observations are the 
black squares linked by the thick black lines, on all three pan- 
els. The 3<x errorbars are shown in the Figure. The model results 
are displayed in colour (blue and red), with the indication of the 
specific shock parameters on each panel, (a) C-type shock mod- 
els: best fit shock parameters «h = 10 4 cirT 3 , no SiO initially 
in the grain mantles; (b) CJ-type shock models: best fit shock 
parameters «h = 10 4 ctrr 3 and u s = 28 km s _1 , with 1% of 
elemental silicon initially in the grain mantles, as SiO; (c) CJ- 
type shock models: best fit shock parameters «h = 10 4 cm 3 and 
y s = 30 km s _1 , with 10% of the silicon initially in the grain 
mantles as SiO. 



carefully. In all cases, as in G08a and G08b, we used an initial 
gas-phase abundance of n(02)/nn = 1.0 x 1 0~ 7 , consistent with 
the 62 observations of the Odin satellite (see Pa gani et al.l l2003: 
lLarsson et alJ 120071) . On the other hand, the initial fractional 
abundance of O2 in chemical equilibrium is 1.0 x 10~ 5 ; we al- 
located the difference between these two values (i.e. most of the 
initial O2) to the grain mantles, consistent with G08a and G08b. 
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Given that the O2 in the mantles is released rapidly into the gas 
phase, by sputtering processes induced by the shock wave, this 
scenario does not differ significantly from that in which all of the 
O2 is assumed to be present initially in the gas phase (the chem- 
ical equilibrium situation). An alternative hypothesis, in which 
the same fraction of the corresponding elemental oxygen is as- 
sumed to be in the form of water ice, rather than O2 ice, was 
considered by G08a and found not to modify significantly the 
results of the models. Indeed, in this case, the rapid erosion of 
the grain mantles during the passage of the shock wave releases 
H2O into the gas phase, and the subsequent reaction of H2O with 
H produces OH, which then reacts with Si to form SiO. A third 
scenario was put to the test in the present study: placing the el- 
emental oxygen in the grain mantles in the form of CO2 ice. 
In this case, CO2 is released quickly into the gas phase during 
the passage of the shock wave, but no chemical reactions con- 
vert it to O2 or OH. When SiO is present in the grain mantles, 
direct sputtering is the dominant source of the gas-phase SiO. 
However, when SiO is absent from the grain mantles, erosion 
of the silicon-bearing grain cores is the mechanism that releases 
silicon into the gas phase, and the deficit in O2 or OH leads to 
lower abundances of gas-phase SiO and weaker SiO line emis- 
sion from low-velocity shocks. In high-velocity shock models, 
or at higher pre-shock densities, O2 and OH are destroyed in the 
gas-phase, through the 02(H,0)OH and OH(H,0)H2 sequence, 
and so the initial composition of the grain mantles is less signif- 
icant. 

If one ignores the observations of H2, it is possible to fit the 
observations of SiO by means of C-type shock models, with- 
out any SiO initially in the grain mantles; cf. G08a. In this case, 
we find that shock models with a pre-shock density of 10 4 cm" 3 
provide better fits than those with pre-shock densities of 10 s or 
10 6 cm" 3 . The typical shock velocities are from 25 to 30 km s _1 , 
with values of the magnetic-field parameter from 0.45 to 2. An 
example of such a fit can be seen in panel (a) of Figure [5] We 
note that, in this case, the assumption that most of the O2 is ini- 
tially in the form of CO2 ice on the grain mantles leads to insuffi- 
cient levels of SiO emission, as compared with the observations. 
However, our aim is to model simultaneously the emission of 
SiO and H2, assuming their emitting regions coincide, to within 
the spatial resolution of the observations. Accordingly, we cal- 
culate the SiO spectrum for the shock models that already fit the 
H2 lines, namely the non-stationary shock models presented in 
Subsection 12.31 Given the the additional flexibility afforded to 
these models by the initial fractional abundance of SiO in the 
grain mantles, we found that most of the CJ-type models that 
already fit the H2 data can reproduce the SiO emission satisfac- 
torily. The results are presented in Figure[5] where panels b) and 
c) correspond to our best fitting models with, respectively, 1% 
and 10% SiO initially in the mantles. We find that no model in 
which SiO is absent from the grain mantles is able to reproduce 
the intensities of the SiO lines simultaneously with those of H2, 
thereby confirming a conclusion of G08b and rendering consid- 
erations of the initial chemical form of oxygen irrelevant. 

We wish to emphasize that a single, plane-parallel CJ-type 
shock wave model is unlikely to provide a completely satisfac- 
tory simulation of the molecular line emission from an object 
as complex as BHR71. In particular, a combination of C-type 
shock models with appropriate parameters and filling factors 
might also be able to reproduce the H2 and SiO line intensities 
and address the problem of fitting the SiO line profiles. The lat- 
ter could be strongly affected by geometrical effects, associated 
with the projection of multiple shock velocities on to the line of 
sight. However, such a study introduces additional free parame- 



ters that are insufficiently constrained by the currently available 
observations. Furthermore, our model does not take into account 
the con sequences of grain-grain collisio ns, whose inclusion (fol- 
lowing iGuillet et al1]2007l I2009L l201lh would have major im- 
plications for computation times. In view of these limitations, 
geometrical and physical, it is not feasible currently to investi- 
gate the nature of the bullet that we observe in the SiO(5^1) line 
profile. 

4. Water emission predictions 

4. 1 . Water modelling 

The abundance of water at each step of the calculation is 
the result mainly of two formation pathways. One is the 
0(H2,H)OH(H2,H)H20 sequence, which takes place in the gas 
phase and becomes signi ficant when the temperature exceeds a 
few hundred kelvin (e.g. ICharnlevll 19971 lAtkinson et alJ l2004). 
The other is the sputtering of water ice that formed in the mantles 
from O and H atoms, which were ads orbed on to the grains in 
the parent dark cloud (see for example (Tie lens & Hagenl fl982). 
In the latter case, the sputtering process is analogous to that con- 
sidered for SiO when present in the grain mantles. 

The modelling of water emission was performed in a sim- 
ilar manner to that of SiO. Again, we used an external ra- 
diative transfer module based on the LVG approximation; 
the molecular data were taken from the LAMDA database 
(http://www.strw.leidenuniv.nl/~moldata/). Regarding the rate 
coefficients for ex citation of H2O by He, we used the data of 
iGreen et aD (Il993l) . For collisions of H2O with H2, we adopted 
the same approac h as FP10 and used the CTMC results of 
iFaure et alJ (120071) . rather than the more precise but less com- 
plete (in terms of their temp e rature r ange) quantal c alcula - 
tion s of iDubernet & Grosieanl d2002l) , iGrosiean etail (|2003), 
and IDubernet et al.1 (120091) . As no data exist for the excitation 
of H2O in collisions with H, we used the same rate coefficients 
as for excitation by ortho-H2, following FP10. 

The more complex structure of H2O, as compared to SiO, 
results in some changes being required to the computational 
method in the radiative transfer module. First, the code must be 
modified to take account of the more complicated radiative se- 
lection rules. Second, in the more complex case of water as com- 
pared to SiO, a value of 10" 10 is required for the convergence 
parameter. Finally, we found it necessary to adapt the expres- 
sion that is used for the escape probability. For a plane-parallel 
geometry, such as is assumed in our shock models, the escape 
probability is 

1 - e" T ^ 

P± = ■ (1) 

T± 

where t ± is the LVG optical depth in the z-direction, normal to 
the shock front. On the other hand, Ne ufeld & Kaufmanl d!993h 
proposed the following approximation to the angle-averaged es- 
cape probability, 



In the case of water, we found that the use of the second for- 
mula gave rise to numerical instabilities and longer computation 
times. Consequently, we adopted the first formula, both when 
solving the equations of statistical equilibrium and when evalu- 
ating the local line emissivity at each point in the shock wave. 
This approach is consistent with that of FP10 and enables us 
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to make direct comparisons with results obtained using their 
method of treating radiative line transfer in shock waves. The 
implications of the different forms of the escape probability for 
the intensities of the rotational transitions of SiO are discussed 
in Appendix A of G08a. We recall that the geometry of the shock 
wave is not well known. 

Before comparing with the method of FP10, we 
benchmarked our LVG code against the RADEX code 
(Ivan der Tak et alJ 12007). at four representative points, cover- 
ing different conditions of temperature, density, and velocity 
gradient, namely (respectively in K, crrT 3 , and s ): (25, 10 4 , 
2 x 10~ 13 ), (50, 10 5 , 10-°), (300, 10 4 , 10~ 12 ), and (467, 
3.71 x 10 5 , 9 x 10 ). For the purpose of this comparison, we 
adopted the same escape probability formula as in the RADEX 
code. For ortho-H20, the computed excitation temperatures 
showed no mean relative discrepancy above 6.3% for the ten 
levels emitting the strongest lines at each of the four points. In 
the case of para-F^O, the relative discrepancy was occasionally 
larger, reaching 43.9%. Although it remains unclear why 
the discrepancies were larger for para- than for ortho-HaO, 
the overall level of agreement with the RADEX code was 
reassuring. 

4.2. Steady-state versus statistical-equilibrium 

In this section, we compare the results of the present calculations 
with those obtained using the approach of FP10. We present the 
intensities of all ortho- and para-transitions, integrated over the 
whole extent of the shock wave. We compared one stationary 
C-shock model, with input parameters «h = 2xl0 4 cm -3 , 
u s - 20 km s _1 , and b = 1, and one non-stationary CJ-shock 
model, with input parameters «h = 10 4 cm~ 3 , o s = 28 km s _1 , 
b = 1.5, and age 500 yrs, which fits both the H2 and SiO ob- 
servations (see Figure [5] panel b). The principal difference be- 
tween the two versions of the shock code is that FP10 included 
the population densities of the levels of water as variables of the 
model, and the differential equations describing the evolution of 
the level populations were solved simultaneously with those for 
the other dynamical and chemical variables, in the steady-state 
approximation. Thus, the water line emissivities were computed 
within the shock code, at each point of the shock profile, and 
could be integrated internally. In the present approach, the level 
population densities (in statistical-equilibrium) and the line in- 
tensities are calculated in a separate LVG module, to which the 
shock profiles are provided as input data. 

In addition to the different approximations used when calcu- 
lating the population densities of the water levels (steady-statfl 
and statistical-equilibriurr0), there are two further, potentially 
significant differences between these two sets of calculations. In 
the steady-state case, the radiative cooling function that is calcu- 
lated at each point includes the contribution of the water emis- 
sion, determined from the level population densities and the line 
escape probabilities. This contribution has an effect on the tem- 
perature profile in regions where H2O is a significant coolant. In 
the statistical-equilibirum assumption, the contribution of wa- 
ter emission to the cooling of t he medium is approximated by 
means of the tables provided bv lNeufeld & Kaufman d!993l) . In 
addition, the effect of adsorption of gas-phase species on the ra- 
dius of the grains is treated differently in the two versions of 
the shock code; this affects the gas-grain interactions in the tail 
of the cooling flow, where adsorption takes place. As a conse- 

2 d/dt = and hence d/df = vd/dz 

3 d/dt = 



quence, the coupling of the gas to the grain temperature, as the 
post-shock region is approached, occurs differently in the two 
versions of the shock code. 

In both the C-type and CJ-type cases, we have three sets of 
LVG calculations. We denote by 'FP10' the results of the steady- 
state (s-s) shock code. Additionally, the dynamical and chemical 
profiles generated by either version of the shock code may be 
combined with a separate, statistical-equilibrium (s-e) treatment 
of the water level populations. The upper panels of Figure|6]com- 
pares the neutral temperature and water fractional abundance 
profiles from the two versions of the shock code. In the case of 
the C-type model, differences are noticeable only in the cooling 
flow, where water is a significant coolant (see top right panel of 
Figure 2). Nonetheless, the differences in the thermal profiles are 
modest, and the water fractional abundances in the two versions 
remain indistinguishable. In the case of the CJ-type model, even 
a small displacement of the location of the J-discontinuity can 
modify the maximum neutral temperature, which in turn has an 
effect on the fractional abundance profiles: see the r.h.s. upper 
panel of Figure [6] 

The lower panels of Figure [6] show the integrated intensities 
of the total of 156 transitions of para-FLO that result from the 
three different calculations. In the case of the C-shock, the sim- 
ilarity of the results is striking and extends to very weak lines. 
The time-scale for the level populations to approach equilibrium 
is given by {n(Yiz)q{J + 1 — > J)]~ l (where q(J +1 — » J) is 
the de-excitation rate coefficient for collisions with molecular 
H2) in the limit of high densities; it is even shorter in the limit 
of low densities, when radiative decay dominates collisional de- 
excitation. To validate the assumption of s-e, one must compare 
this population transfer time-scale with a dynamical time-scale, 
such as the flow time of the neutral particles through the shock 
region, that corresponds to the width of the neutral temperature 
profile shown in the upper panels. Adopting the representative 
values of 10 4 ctrT 3 for the density and 10~ 10 crrT 3 for the colli- 
sional de-excitation rate coefficient results in a population trans- 
fer time-scale of 0.03 yr, much less than the flow time of the neu- 
tral particles through the C-type shock considered. In the case of 
the reference CJ-type shock model, this condition is verified for 
most of the transitions, as the overall agreement between light 
and dark blue open circles in Figure [6] confirms. For the bright- 
est transitions, our independent LVG calculation is consistent 
with that denoted 'FP10'. In this case, the discrepancies with 
the red filled circles are larger, owing to the differences in the 
water abundance profiles, which are attributable in turn to the 
temperature differences arising from: the small displacement of 
J-discontinuity; the treatment of the cooling due to water; and 
the strength of the grain-gas interaction as the post-shock region 
is approached. We note that the pattern created by each set of 
points is similar for the brightest, and hence observable, transi- 
tions, and that the pattern is shifted vertically from one set of 
results to the other. Such a shift is indistinguishable from that 
arising from the uncertainty in the filling factor, when compar- 
isons are made with observations. 

We conclude from the above analysis that the use of the s-e 
approximation is justified in the present context of the study of 
water emission. 

4.3. Water emission predictions for BHR71 

As a preamble to this section - aimed at providing predictions 
of water line intensities for use in the interpretation of obser- 
vations - we study the evolution of the optical depth of the 
fundamental transitions of both 0-H2O and p-FLO in each of 
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Fig. 6. Upper panels: comparison between profiles obtained when the computation of the water line emission is internal to the shock 
code (steady-state approximation: blue curves) or external to the shock code (statistical equilibrium approximation: red curves); the 
temperature of the neutral fluid (continuous line) and the fractional abundance of water (dashed lines) are shown for our two ref- 
erence, C-type (left-hand panel), and CJ-type (right-hand panel) shock models. In both cases, the shock parameters are displayed 
in the corresponding panels. Lower panels: integrated line intensities for all para-water transitions; steady-state approximation to 
the computation of the water emission (dark blue open circles); same shock model, but with the statistical-equilibrium approxima- 
tion to the computation of the water emission (large light-blue open circles); present shock model, with the statistical-equilibrium 
assumption (small red filled circles). 



our reference models. This evolution, together with the total op- 
tical depths, integrated through the shock wave, are shown in 
Figure Q Owing to the one-dimensional nature of our approach, 
we probably overestimate the optical depths and find that the 
transitions are optically thick in most parts of our shock mod- 
els; this is especially the case of the C-type shock, as the water 
accumulates in the large post-shock region. The example of the 
CJ-type shock shows the spread of values of the optical depths 
in the emitting region - low when water is dissociated in the 
high-temperature regime, followed by higher values when water 
reforms. This figure underlines the requirement for an adequate 
treatment of the radiative transfer. 

As already mentioned, H2, SiO and H2O are believed to have 
coincident regions of emission, at least to the resolu tion of cur- 
rently available observations (e.g. Nisi ni et al.ll2.Q10r) . For com- 
pleteness, we provide predictions of water line intensities for all 
the models that were found to fit the H2 pure rotational excitation 



diagram acceptably at position of the SiO knot. Given the rela- 
tive uniformity of the shock conditions in the inner parts of the 
outflow, revealed by the similarity of the H2 excitation diagrams 
at the different positions considered, our models should be ap- 
plicable to water observations at any position. Stationary, C-type 
shock models are excluded, for the reasons given in Section [33l 

The results are given in Tables lBTI and lTJ!2l for o-FbO and p- 
H2O, respectively. The transitions that appear in the Tables are 
the ones that will be the target of Herschel observations. The first 
three transitions in each table will be observed by HIFI, and the 
remainder (nine transitions of 0-H2O and three of P-H2O) will be 
observed by PACS, in the framework of the WISH key-program 
(van Dishoeck et al The calculations comprise the entire 

set of transitions that are listed on the LAMDA website, and the 
complete results are available upon request. 

The first general comment concerns the dependence of the 
water emission line intensities on the pre-shock density. We find 
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Fig. 7. Evolution of the computed local optical depths of two transitions as a function of the flow time of the neutral fluid in our 
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that the integrated intensities calculated for a pre-shock density 
of 10 5 cm 4 are systematically around one order of magnitude 
larger than those computed in models with pre-shock densities 
of 10 4 crrT 3 , regardless of the values of the other parameters, 
such as the shock velocity or magnetic-field strength. In addi- 
tion, and with all other parameters remaining equal, we find that 
the computed integrated intensities vary by factors of 2-3, and 
sometimes more, over the range of shock ages considered; this 
variation may be sufficient, in certain circumstances, to enable 
the observations to constrain the age of the shock wave. 

Finally, we consider the computed data pertaining to the 
eight shock models that provide the best simultaneous fits for H2 
and SiO; these models are shown in Figure [5] in panels (b) and 
(c). The results of the calculations of the integrated intensities 
of the brightest o- and P-H2O transitions, likely to be targetted 
by HIFI and PACS, are shown in Figure [8] The filling factor is 
assumed to be the same and equal to 1 for all the transitions and 
too small to discriminate between these models. Nevertheless, 
forthcoming Herschel observations should enable us to enhance 
our understanding of the chemistry of water in interstellar shock 
waves. 

Water is believed to be the main interstellar reservoir of oxy- 
gen, which is the most abundant element in the universe after H 
and He and plays a major role in interstellar chemistry. Given 
its rotational energy level structure, water is an important gas 
coolant and a convenient tracer of shocked regions. Forthcoming 
observations, through their comparison with the predictions of 
shock models, should enable us to investigate both the chemical 
and the physical processes at work in bipolar outflows. 



5. Concluding remarks 

We have presented a detailed analysis of observations at a few 
positions in the the bipolar outflow BHR71, assuming that the 
emission from these regions arises in shock waves. In two cases, 
namely the Herbig-Haro objects HH320 and HH321, we are able 
to make use of both the recent Spitzer observations of pure ro- 
tational transitions and previous, ground-based observations of 
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and P-H2O transitions are indicated by circles and squares, re- 
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density of 10 4 cm 4 , a magnetic-field parameter of 1.5, a shock 
velocity of 28 km s , and an age of 500 years. 
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rovibrational transitions of molecular hydrogen. In order to com- 
pare with the predictions of our shock models, we used excita- 
tion diagrams. The physical conditions were found to be similar 
in both cases, with a pre-shock density of 10 4 crrT 3 and moder- 
ate shock velocities, in the range 15-35 km s _1 . The only shock 
models that fit these data satisfactorily are non-stationary, with 
ages that depend on the magnetic-field strength and do not ex- 
ceed 2000 years. Three other positions in the upper and lower 
lobes of the outflow were studied, using only the pure rotational 
transitions of H2 observed by Spitzer. Our analysis suggests that 
the excitation conditions at these positions are similar to those 
in the HH objects and can only be reproduced by non-stationary 
models. 

We present a new map of the SiO(5^) transition in the inner 
part of the bipolar outflow BHR71, together with new observa- 
tions of the SiO(6-5) and SiO(8-7) lines, at one position that 
is also covered by the Spitzer observations and at another posi- 
tion that lies outside of the Spitzer maps. At the first position, 
referred to as 'the SiO knot', we make use of the combination 
of data (pure rotational transitions of both H2 and SiO) in order 
to constrain the parameters of our shock models more tightly 
than is possible with the rotational transitions of H2 alone. The 
modelling based on H2 leads to the conclusion as for the other 
positions in the outflow, namely non-stationary shocks provide 
the best fits. In particular, shocks with a pre-shock density of 
10 5 crrT 3 , low velocities (10-20 km s _I ) and ages (less than 
1000 years) were found to be good candidates. Including the SiO 
emission places constraints on the initial repartition of silicon- 
bearing material but is of little help in discriminating between 
the models which fit the H2 rotational excitation diagram. 

We have validated our approach to calculating the water line 
intensities, which uses the results of shock model as input to 
a separate radiative transfer code, based on the LVG approx- 
imation. To this end, we provide a full comparison between 
the statistical-equilibrium and the steady-state approximations, 
based on reference shock models of the stationary, C-type, and 
non-stationary, CJ-type. We make predictions of water emission 
line intensities for the SiO knot in the BHR71 outflow. As we 
expect emission by water to coincide with that by SiO and H2, 
we provide the results of the LVG calculations of the water line 
intensities for models that fit satisfactorily both the SiO and H2 
emission^ Given that the excitation conditions are relatively uni- 
form in the internal parts of the outflow, we consider that our 
calculations can be used to compare with observations of water 
throughout this region. 

We shall discuss in a future publication ground-based ob- 
servations of CO transitions, undertaken with the APEX facil- 
ity, that could bring further insights into the physical conditions 
in the BHR7 1 outflow. Complementary observations of sulphur- 
bearing molecules and other shock tracers, such as HNCO, are 
planned. One of the aims of these studies is to investigate the 
suitability of 'chemically active outflows', such as BHR71, for 
observations with the new generation of instruments, such as 
ALMA. 
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4 On request, we can provide water line intensities for all the shock 
models that fit the H2 observations. 
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Appendix A: Model constraints for the HH objects 
and the SiO knot 
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Table A.l. Parameters of the models that fit the H2 observations 
of HH320AB and HH321AB (separated by the double horizon- 
tal line). The shock velocities are in km s _1 and the ages in years. 
Ages between agel and age2 provide acceptable fits to the ob- 
servations of Ht. 
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Table A.2. Parameters of the models that fit the H2 observations 
of the SiO knot. The shock velocities are in km s _1 and the ages 
in years. Ages between agel and age2 provide acceptable fits to 
the observations of H2. 
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0.60 


22 


275 


365 


CJ 


4 


0.75 


15 


710 


1120 


CJ 


4 


0.75 


18 


375 


695 


CJ 


4 


0.75 


20 


365 


530 


CJ 


4 


0.75 


22 


225 


450 


CJ 


4 


1.00 


20 


490 


740 


CJ 


4 


1.00 


22 


370 


465 


CJ 


4 


1.00 


25 


290 


345 


CJ 


4 


1.25 


22 


585 


670 


CJ 


4 


1.25 


25 


425 


530 


CJ 


4 


1.25 


28 


325 


465 


CJ 


4 


1.50 


22 


830 


1005 


CJ 


4 


1.50 


25 


605 


755 


CJ 


4 


1.50 


28 


500 


625 


CJ 


4 


1.50 


30 


405 


825 


CJ 


4 


1.75 


22 


1115 


1260 


CJ 


4 


1.75 


25 


855 


920 


CJ 


4 


1.75 


28 


660 


890 


CJ 


4 


1.75 


30 


570 


670 


CJ 


4 


2.00 


28 


835 


920 


CJ 


4 


2.00 


30 


730 


775 


CJ 


5 


0.30 


10 


16 


67 


CJ 


5 


0.30 


12 


23 


69 


CJ 


5 


0.45 


10 


51 


115 


CJ 


5 


0.45 


12 


58 


76 


CJ 


5 


0.45 


15 


24 


55 


CJ 


5 


0.45 


18 


13 


35 


CJ 


5 


0.60 


10 


52 


155 


CJ 


5 


0.60 


12 


77 


120 


CJ 


5 


0.60 


15 


53 


72 


CJ 


5 


0.60 


20 


26 


62 


CJ 


5 


0.75 


10 


115 


250 


CJ 


5 


0.75 


12 


100 


165 


CJ 


5 


0.75 


15 


86 


104 


CJ 


5 


1.00 


10 


385 


410 


CJ 


5 


1.00 


12 


205 


255 


CJ 


5 


1.00 


15 


160 


175 


CJ 


5 


1.25 


15 


205 


265 


CJ 


5 


1.50 


10 


525 


775 


CJ 


5 


1.50 


15 


335 


370 


CJ 


5 


1.75 


10 


505 


915 


CJ 


5 


1.75 


15 


315 


465 


CJ 


5 


1.75 


18 


305 


340 


CJ 


5 


2.00 


18 


435 


440 
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Table B.l. Predicted integrated intensities of the 0-H2O lines to be targeted by the HIFI and PACS receivers of the Herschel telescope. The units are K km s _1 . The models 
considered are those which best fit the H2 pure rotational excitation diagram at the position of the SiO knot. A model referred to as 'n4b0.45vl5t335' has the following set of 
input parameters: a pre-shock density of 10 4 cirT 3 , a magnetic-field parameter of 0.45, a shock velocity of 15 km s _1 , and an age of 335 years. 



transition name 


A 10 " iOl 


-x , _ 1 n „ 
J12 ~ JQ3 




^21 " ^12 


7 .„ 1 ,-, , 


J03 " ^12 


■^30 " J21 


4l4 " 3()3 


^21 " A 10 


7()7 - 616 


8 . „ 7„ 

°18 " '07 


?Ql - 018 


frequency (GHz) 


JJO.VjO 




1 1 3 J. IZ 1 


lOOl.UUo 


1 AAQ on^ 
looy.yuj 


1 71 A 77H 
1 / 10. / /U 


9 1 OA 1 




0771 Q77 


41 AA 8^7 




^77A 

JZ /O. JZU 


n4b0.45vl5t335 


5.0E+00 


3.6E+00 


4.7E-01 


6.4E-01 


2.1E+00 


1.7E+00 


1.IE+00 


3.9E-01 


2.8E-02 


9.6E-03 


4.0E-03 


1.5E-03 


n4b0.45vl5t555 


6.6E+00 


4.5E+00 


5.9E-01 


7.8E-01 


2.7E+00 


2.2E+00 


1.4E+00 


5.3E-01 


3.9E-02 


1.5E-02 


6.2E-03 


2.4E-03 


n4b0.45vl5t645 


7.3E+00 


4.8E+00 


6.3E-01 


8.3E-01 


2.9E+00 


2.3E+00 


1.5E+00 


5.8E-01 


4.2E-02 


1.7E-02 


7.1E-03 


2.7E-03 


n4b0.45vl8t350 


7.6E+00 


5.8E+00 


1.0E+00 


1.2E+00 


3.8E+00 


4.4E+00 


2.0E+00 


9.7E-01 


1.0E-01 


2.3E-02 


9.2E-03 


3.4E-03 


n4b0.45vl8t385 


7.9E+00 


5.9E+00 


1.1E+00 


1.3E+00 


3.9E+00 


4.5E+00 


2.1E+00 


1.0E+00 


1.0E-01 


2.4E-02 


9.8E-03 


3.6E-03 


n4b0.45vl8t425 


8.2E+00 


6.1E+00 


1.1E+00 


1.3E+00 


4.0E+00 


4.6E+00 


2.2E+00 


1.0E+00 


1.1E-01 


2.6E-02 


1.1E-02 


3.9E-03 


n4b0.45v20t210 


7.6E+00 


6.3E+00 


1.3E+00 


1.5E+00 


4.4E+00 


5.9E+00 


2.4E+00 


1.3E+00 


1.5E-01 


2.9E-02 


1.1E-02 


3.9E-03 


n4b0.45v20t315 


8.8E+00 


7.0E+00 


1.5E+00 


1.6E+00 


4.8E+00 


6.4E+00 


2.7E+00 


1.4E+00 


1.6E-01 


3.5E-02 


1.3E-02 


4.9E-03 


n4b().45v20t600 


1.1E+01 


7.5E+00 


1.4E+00 


1.6E+00 


5.0E+00 


6.0E+00 


2.8E+00 


1.4E+00 


1.5E-01 


3.9E-02 


1.6E-02 


6.0E-03 



n4b0.60vl5t640 


5.9E+00 


3.8E+00 


4.7E-01 


6.5E-01 


2.2E+00 


1.7E+00 


1 


.1E+00 


4.0E-01 


2.9E- 


■02 


LIE- 


■02 


4.6E- 


■03 


1.7E- 


■03 


n4b0.60vl5t655 


6.0E+00 


3.9E+00 


4.8E-01 


6.6E-01 


2.3E+00 


1.8E+00 


1 


.2E+00 


4.1E-01 


3.0E- 


■02 


LIE- 


■02 


4.8E- 


■03 


1.8E- 


■03 


n4b0.60vl5t885 


7.7E+00 


4.6E+00 


5.5E-01 


7.5E-01 


2.6E+00 


2.0E+00 


1 


.4E+00 


5.1E-01 


3.6E- 


■02 


1.5E- 


■02 


6.4E- 


■03 


2.4E- 


■03 


n4b0.60v 181370 


6.6E+00 


4.9E+00 


8.8E-01 


1.1E+00 


3.2E+00 


3.6E+00 


1 


.7E+00 


7.6E-01 


8.1E- 


■02 


1.7E- 


■02 


6.7E- 


■03 


2.4E- 


■03 


n4b0.60vl8t515 


8.1E+00 


5.7E+00 


9.7E-01 


1.2E+00 


3.7E+00 


4.0E+00 


2 


.OE+00 


8.9E-01 


9.1E- 


■02 


2.2E- 


■02 


9.0E- 


■03 


3.4E- 


■03 


n4b0.60vl8t730 


1.0E+01 


6.7E+00 


1.1E+00 


1.3E+00 


4.2E+00 


4.5E+00 


2 


.3E+00 


1.1E+00 


1.0E- 


-01 


3.0E- 


■02 


1.2E- 


■02 


4.7E- 


■03 


n4b0.60v20t260 


7.1E+00 


5.7E+00 


1.2E+00 


1.4E+00 


3.9E+00 


5.3E+00 


2 


.2E+00 


1.1E+00 


1.3E- 


■01 


2.3E- 


■02 


8.6E- 


■03 


3.1E- 


■03 


n4b().60v20t330 


8.0E+00 


6.3E+00 


1.3E+00 


1.5E+00 


4.3E+00 


5.7E+00 


2 


.4E+00 


1.2E+00 


1.4E- 


■01 


2.7E- 


■02 


1.0E- 


■02 


3.8E- 


■03 


n4b0.60v20t560 


1.0E+01 


7.3E+00 


1.3E+00 


1.6E+00 


4.8E+00 


5.8E+00 


2 


.7E+00 


1.3E+00 


1.4E- 


-01 


3.5E- 


■02 


1.4E- 


■02 


5.3E- 


■03 


n4b0.60v22tl95 


8.5E+00 


7.7E+00 


1.9E+00 


2.0E+00 


5.4E+00 


8.8E+OO 


3 


.2E+00 


1.9E+00 


2.4E- 


■01 


4.5E- 


■02 


1.5E- 


■02 


5.3E- 


■03 


n4b0.60v22t275 


9.9E+00 


8.4E+00 


2.0E+00 


2.1E+00 


5.9E+00 


9.4E+00 


3 


.6E+00 


2.0E+00 


2.6E- 


■01 


4.9E- 


■02 


1.7E- 


■02 


6. IE 


■03 


n4b().60v22t365 


1.1E+01 


8.6E+00 


2.0E+00 


2.1E+00 


5.9E+00 


8.9E+00 


3 


.6E+00 


1.9E+00 


2.4E- 


■01 


4.8E- 


■02 


1.8E- 


■02 


6.5E 


■03 



n4b0.75vl5t710 


4.8E+00 


3.1E+00 


3.5E-01 


5.1E-01 


1 


7E+00 


1 


3E+00 


8.5E-01 


2.9E-01 


2. 


0E- 


■02 


7.2E- 


-03 


3. 


IE 


-03 


1 


IE 


-03 


n4b0.75vl5t830 


5.7E+00 


3.5E+00 


4.0E-01 


5.6E-01 


1 


9E+00 


1. 


4E+00 


9.7E-01 


3.3E-01 


2. 


3E- 


■02 


9.0E- 


■03 


3. 


8E- 


■03 


1 


4E- 


-03 


n4b0.75vl5tll20 


7.4E+00 


4.0E+00 


4.2E-01 


6.0E-01 


2. 


2E+00 


1 


5E+00 


1.1E+00 


3.9E-01 


2. 


4E- 


•02 


LIE- 


■02 


5. 


0E- 


-03 


1 


9E- 


■03 


n4b0.75vl8t375 


4.8E+00 


3.7E+00 


6.2E-01 


7.9E-01 


2. 


3E+00 


2. 


5E+00 


1.2E+00 


4.9E-01 


5. 


3E- 


■02 


9.5E- 


■03 


3. 


7E- 


■03 


1 


3E 


-03 


n4b0.75vl8t510 


6.4E+00 


4.5E+00 


7.7E-01 


9.6E-01 


2. 


9E+00 


3. 


1E+00 


1.5E+00 


6.4E-01 


6 


8E- 


■02 


1.4E- 


■02 


5. 


8E- 


■03 


2. 


IE 


-03 


n4b0.75vl8t695 


8.5E+00 


5.6E+00 


8.8E-01 


1.1E+00 


3. 


5E+00 


3. 


6E+00 


1.8E+00 


8.0E-01 


7. 


8E- 


■02 


2.1E- 


■02 


8. 


5E- 


-03 


3. 


2E- 


-03 


n4b0.75v20t365 


7.1E+00 


5.5E+00 


1.1E+00 


1.3E+00 


3. 


7E+00 


4. 


9E+00 


2.1E+00 


9.9E-01 


1 


2E- 


-01 


2.1E- 


■02 


7. 


8E- 


-03 


2. 


8E- 


-03 


n4b0.75v20t430 


7.9E+00 


5.9E+00 


1.2E+00 


1.4E+00 


3. 


9E+00 


5. 


1E+00 


2.2E+00 


1.1E+00 


1 


3E- 


-01 


2.4E- 


■02 


9. 


1E 


■03 


3. 


3E 


■03 


n4b0.75v20t530 


9.1E+00 


6.5E+00 


1.2E+00 


1.5E+00 


4. 


3E+00 


5, 


4E+00 


2.4E+00 


1.2E+00 


1 


3E- 


-01 


2.8E- 


■02 


1 


IE 


02 


4. 


IE 


■03 


n4b0.75v22t225 


7.1E+00 


6.4E+00 


1.5E+00 


1.6E+00 


4. 


5E+00 


7. 


2E+00 


2.6E+00 


1.5E+00 


1 


9E 


01 


3.3E- 


-02 


1 


IE 


02 


3. 


7E- 


-03 


n4b0.75v22t285 


8.4E+00 


7.3E+00 


1.8E+00 


1.9E+00 


5. 


0E+00 


8 


2E+00 


3.1E+00 


1.7E+00 


2, 


2E- 


01 


3.8E- 


-02 


1 


3E- 


■02 


4. 


5E- 


-03 


n4b0.75v22t450 


1.1E+01 


8.2E+00 


1.8E+00 


2.0E+00 


5. 


5E+00 


8. 


1E+00 


3.4E+00 


1.7E+00 


2. 


2E- 


-01 


4.3E- 


■02 


1 


6E 


■02 


5. 


8E- 


-03 



Table B.l. (continued) Predicted integrated intensities of the 0-H2O lines to be targeted by the HIFI and PACS receivers of the Herschel telescope. The units are K km s _1 . The 
models considered are those which best fit the H2 pure rotational excitation diagram at the position of the SiO knot. A model referred to as 'n4b0.45vl5t335' has the following 
set of input parameters: a pre-shock density of 10 4 cm 3 , a magnetic-field parameter of 0.45, a shock velocity of 15 km s _1 , and an age of 335 years. 



transition name 


A 10 - hn 


■3 , _ -2„„ 
J 12 ~ JQ3 


J12 ~ *21 


, 7,„ 
^21 " ^12 


7,„ 1 , 




■^30 " J21 


4l4 - 3o3 


7„ , 1 , „ 
^21 " 110 


7o7 - 616 


C . „ 7„ 

0I8 " '07 


yoi ■ s i8 


frequency (GHz) 


jjo.Vjd 




1 1 3 J. 1Z 1 


lOOl.UUo 


1 AAQ on^ 
looy.yuj 


1 71 A 77fl 
1 1 10. / /U 


9 1 OA 1 




7771 Q77 

Z / / j.y 1 1 


41 £A 8^7 




^77A S7fl 
JZ /O. JZU 


n4bl.00v20t490 


5.3E+00 


4.1E+00 


8.IE-01 


9.7E-01 


2.7E+00 


3.5E+00 


1.5E+00 


6.5E-01 


8.0E-02 


1.2E-02 


4.4E-03 


1.5E-03 


n4bl.00v20t605 


7.4E+00 


5.3E+00 


9.9E-01 


1.2E+00 


3.4E+00 


4.2E+00 


1.9E+00 


8.5E-01 


l.OE-01 


1.9E-02 


7.3E-03 


2.6E-03 


n4bl.00v20t740 


8.9E+00 


6.0E+00 


1.1E+00 


1.3E+00 


3.8E+00 


4.5E+00 


2.1E+00 


9.7E-01 


1.1E-01 


2.3E-02 


9.3E-03 


3.4E-03 


n4bl.00v22t370 


6.0E+00 


5.3E+00 


1.2E+00 


1.3E+00 


3.6E+00 


5.6E+00 


2.1E+00 


1.1E+00 


1.5E-01 


2.1E-02 


7.0E-03 


2.3E-03 


n4bl.00v22t435 


7.5E+00 


6.2E+00 


1.5E+00 


1.6E+00 


4.2E+00 


6.6E+00 


2.6E+00 


1.3E+00 


1.8E-01 


2.6E-02 


9.2E-03 


3.2E-03 


n4bl.00v22t465 


8.0E+00 


6.4E+00 


1.5E+00 


1.6E+00 


4.3E+00 


6.7E+00 


2.6E+00 


1.3E+00 


1.8E-01 


2.8E-02 


9.8E-03 


3.4E-03 


n4bl.00v25t290 


7.3E+00 


6.8E+00 


1.7E+00 


1.8E+00 


4.9E+00 


8.6E+00 


2.9E+00 


1.8E+00 


2.3E-01 


5.2E-02 


1.7E-02 


5.3E-03 


n4bl.00v25t320 


9.1E+00 


8.6E+00 


2.4E+00 


2.3E+00 


6.0E+00 


1.2E+01 


3.9E+00 


2.3E+00 


3.3E-01 


7.3E-02 


2.2E-02 


7.1E-03 


n4bl.00v25t345 


1.1E+01 


9.9E+00 


2.8E+00 


2.7E+00 


6.8E+00 


1.4E+01 


4.6E+00 


2.7E+00 


4.1E-01 


8.6E-02 


2.6E-02 


8.3E-03 



n4bl.25v22t585 


5.7E+00 


4.7E+00 


1.0E+00 


1.2E+00 


3.1E+00 


4.6E+00 


1.8E+00 


8.6E-01 


1.2E-01 


1.5E-02 


5.2E-03 


1.7E-03 


n4bl.25v22t625 


6.7E+00 


5.2E+00 


1.2E+00 


1.3E+00 


3.5E+00 


5.1E+00 


2.1E+00 


9.7E-01 


1.3E-01 


1.8E-02 


6.5E-03 


2.3E-03 


n4bl.25v22t670 


7.3E+00 


5.6E+00 


1.2E+00 


1.4E+00 


3.7E+00 


5.3E+00 


2.2E+00 


1.0E+00 


1.4E-01 


2.0E-02 


7.4E-03 


2.6E-03 


n4bl.25v25t425 


7.5E+00 


7.3E+00 


2.0E+00 


2.0E+00 


5.1E+00 


1.0E+01 


3.3E+00 


1.9E+00 


2.8E-01 


6.1E-02 


1.8E-02 


5.6E-03 


n4bl.25v25t445 


8.2E+00 


7.9E+00 


2.2E+00 


2.2E+00 


5.4E+00 


1.1E+01 


3.6E+00 


2.1E+00 


3.2E-01 


6.6E-02 


2.0E-02 


6.0E-03 


n4bl.25v25t530 


1.1E+01 


1.0E+01 


2.8E+00 


2.8E+00 


6.7E+00 


1.3E+01 


4.7E+00 


2.6E+00 


4.5E-01 


7.6E-02 


2.3E-02 


7.3E-03 


n4bl.25v28t325 


7.0E-01 


2.0E-01 


7.4E-03 


6.2E-03 


9.1E-02 


6.9E-03 


4.8E-02 


2.2E-02 


6.0E-05 


1.3E-03 


6.4E-04 


2.7E-04 


n4bl.25v28t465 


4.8E+00 


2.0E+00 


5.9E-02 


1.5E-01 


9.3E-01 


7.9E-02 


5.3E-01 


2.0E-01 


8.1E-04 


1.2E-02 


5.9E-03 


2.5E-03 



n4bl.50v22t830 


4.9E+00 


4.0E+00 


8.2E-01 


9.6E-01 


2.6E+00 


3.6E+00 


1.4E+00 


6.3E-01 


8.5E-02 


1.0E-02 


3.7E-03 


1.3E-03 


n4bl.50v22t930 


6.9E+00 


5.0E+00 


1.0E+00 


1.2E+00 


3.2E+00 


4.3E+00 


1.9E+00 


8.2E-01 


1.1E-01 


1.6E-02 


6.0E-03 


2.1E-03 


n4bl.50v22tl005 


8.2E+00 


5.6E+00 


1.1E+00 


1.3E+00 


3.5E+00 


4.6E+00 


2.1E+00 


9.0E-01 


1.1E-01 


2.0E-02 


7.5E-03 


2.7E-03 


n4bl.50v25t605 


7.1E+00 


6.9E+00 


1.9E+00 


1.9E+00 


4.7E+00 


9.3E+00 


3.1E+00 


1.7E+00 


2.6E-01 


4.6E-02 


1.4E-02 


4.1E-03 


n4bl.50v25t655 


8.4E+00 


7.7E+00 


2.1E+00 


2.1E+00 


5.2E+00 


1.0E+01 


3.5E+00 


1.9E+00 


3.1E-01 


5.0E-02 


1.5E-02 


4.6E-03 


n4bl.50v25t755 


1.1E+01 


8.9E+00 


2.3E+00 


2.4E+00 


5.9E+00 


1.1E+01 


4.0E+00 


2.1E+00 


3.5E-01 


5.1E-02 


1.7E-02 


5.5E-03 


n4bl.50v28t500 


2.4E+00 


1.0E+00 


4.1E-02 


7.8E-02 


4.6E-01 


5.2E-02 


2.7E-01 


8.4E-02 


4.7E-04 


3.9E-03 


1.9E-03 


7.5E-04 


n4bl.50v28t520 


2.6E+00 


1.1E+00 


4.0E-02 


7.9E-02 


4.8E-01 


4.8E-02 


2.8E-01 


9.0E-02 


4.5E-04 


4.5E-03 


2.1E-03 


8.6E-04 


n4bl.50v28t625 


4.7E+00 


2.0E+00 


5.6E-02 


1.5E-01 


8.9E-01 


7.7E-02 


5.2E-01 


1.9E-01 


7.7E-04 


1.1E-02 


5.3E-03 


2.2E-03 


n4bl.50v30t405 


1.8E+00 


6.6E-01 


2.4E-02 


4.2E-02 


2.9E-01 


2.5E-02 


1.6E-01 


5.7E-02 


2.6E-04 


2.8E-03 


1.4E-03 


5.6E-04 


n4bl.50v30t435 


2.2E+00 


8.0E-01 


2.6E-02 


5.0E-02 


3.5E-01 


2.9E-02 


2.0E-01 


7.5E-02 


2.9E-04 


4.2E-03 


2.1E-03 


8.7E-04 


n4bl.50v30t825 


1.3E+01 


6.0E+00 


1.6E-01 


4.8E-01 


2.8E+00 


2.5E-01 


1.7E+00 


6.3E-01 


2.6E-03 


3.9E-02 


1.9E-02 


8.0E-03 



Table B.l. (continued) Predicted integrated intensities of the 0-H2O lines to be targeted by the HIFI and PACS receivers of the Herschel telescope. The units are K km s _1 . The 
models considered are those which best fit the H2 pure rotational excitation diagram at the position of the SiO knot. A model referred to as 'n4b0.45vl5t335' has the following 
set of input parameters: a pre-shock density of 10 4 cm 3 , a magnetic-field parameter of 0.45, a shock velocity of 15 km s _1 , and an age of 335 years. 



transition name 
frequency (GHz) 


1 .r. 1„, 

r 10 - 101 

JJO.VJD 


312 " 303 
luy/.Joj 


7~, 

J12 - Z21 

1 1 C-3 1 77 
1 1 J J. 1Z / 


7~, 7.~ 
Z21 - Z12 

lODl.UUo 


7,~ 

Z12 " iOl 
1 AAQ QO^ 

looy.yuj 


7,~ 

J03 " Z12 
1 71 £ 77H 

1 / 10. / /u 


J30 " J21 

7 1 OA 1 A& 
Z1V0.J40 


4l4 - 3o3 

7AAH A1A 
ZD4U.4 /4 


7,, 1 
Z2i - r 10 

7771 Q77 
Z/ Ij.yl 1 


7o7 - 616 

A\ AA 8^7 
4100.0JZ 


C . „ 7- 
018 " '07 

ATiA 7QA 


yoi - 018 

C77A ^7n 

JZ /O. JZU 


n4bl.75v22tlll5 


4.5E+00 


3.5E+O0 


7.0E-01 


8.3E-01 


2.2E+00 


3.0E+00 


1.2E+00 


5.2E-01 


7.0E-02 


8.5E-03 


3.0E-03 


1.0E-03 


n4bl.75v22tll80 


5.5E+00 


4.0E+00 


7.6E-01 


9.3E-01 


2.5E+00 


3.2E+00 


1.4E+00 


5.8E-01 


7.5E-02 


1.1E-02 


4.0E-03 


1.4E-03 


n4bl.75v22tl260 


6.8E+00 


4.6E+00 


8.5E-01 


1.0E+00 


2.9E+00 


3.6E+00 


1.6E+00 


6.8E-01 


8.4E-02 


1.4E-02 


5.3E-03 


1.8E-03 


n4bl.75v25t855 


7.3E+00 


6.5E+00 


1.7E+00 


1.8E+00 


4.4E+00 


8.3E+00 


2.9E+00 


1.5E+00 


2.4E-01 


3.3E-02 


1.0E-02 


3.2E-03 


n4bl.75v25t885 


7.9E+00 


6.8E+00 


1.8E+00 


1.9E+00 


4.6E+00 


8.4E+00 


3.0E+00 


1.6E+00 


2.5E-01 


3.4E-02 


1.0E-02 


3.4E-03 


n4bl.75v25t920 


8.7E+00 


7.2E+00 


1.8E+00 


1.9E+00 


4.8E+00 


8.5E+00 


3.2E+00 


1.6E+00 


2.6E-01 


3.5E-02 


1.1E-02 


3.7E-03 


n4bl.75v28t660 


2.0E+00 


9.0E-01 


4.3E-02 


7.8E-02 


4.1E-01 


6.7E-02 


2.5E-01 


6.7E-02 


5.4E-04 


2.5E-03 


1.1E-03 


4.3E-04 


n4bl.75v28t705 


2.9E+00 


1.2E+00 


4.4E-02 


9.2E-02 


5.3E-01 


5.8E-02 


3.1E-01 


9.8E-02 


5.1E-04 


4.7E-03 


2.2E-03 


8.7E-04 


n4bl.75v28t890 


1.5E+01 


1.2E+01 


2.7E+00 


2.9E+00 


7.4E+00 


1.3E+01 


5.1E+00 


2.7E+00 


4.2E-01 


7.7E-02 


2.5E-02 


8.5E-03 


n4bl.75v30t570 


1.9E+00 


7.1E-01 


2.5E-02 


4.7E-02 


3.0E-01 


2.7E-02 


1.7E-01 


5.8E-02 


2.6E-04 


2.9E-03 


1.4E-03 


5.5E-04 


n4bl.75v30t600 


2.6E+00 


9.8E-01 


3.0E-02 


6.6E-02 


4.3E-01 


3.5E-02 


2.4E-01 


8.9E-02 


3.5E-04 


4.9E-03 


2.4E-03 


9.7E-04 


n4bl.75v30t670 


4.3E+00 


1.8E+00 


5.0E-02 


1.3E-01 


8.1E-01 


7.0E-02 


4.7E-01 


1.8E-01 


7.0E-04 


1.0E-02 


5.0E-03 


2.0E-03 



n4b2.00v28t835 


1.8E+00 


1.1E+00 


8.3E-02 


1.3E-01 


5.5E-01 


2.3E-01 


3.1E-01 


8.3E-02 


2.3E-03 


2.1E-03 


8.4E-04 


2.8E-04 


n4b2.00v28t880 


2.8E+00 


1.5E+00 


1.1E-01 


1.8E-01 


7.7E-01 


3.1E-01 


4.2E-01 


1.3E-01 


3.4E-03 


4.0E-03 


1.7E-03 


6.3E-04 


n4b2.00v28t920 


5.2E+00 


3.3E+00 


3.9E-01 


5.5E-01 


1.9E+00 


1.5E+00 


9.5E-01 


3.8E-01 


3.2E-02 


1.0E-02 


4.2E-03 


1.5E-03 


n4b2.00v30t730 


1.7E+00 


6.3E-01 


2.3E-02 


4.4E-02 


2.6E-01 


2.5E-02 


1.5E-01 


4.6E-02 


2.4E-04 


2.0E-03 


9.2E-04 


3.5E-04 


n4b2.00v30t755 


2.2E+00 


8.4E-01 


3.0E-02 


6.1E-02 


3.6E-01 


3.5E-02 


2.1E-01 


6.8E-02 


3.2E-04 


3.2E-03 


1.5E-03 


5.9E-04 


n4b2.00v30t775 


2.6E+00 


9.9E-01 


3.1E-02 


7.0E-02 


4.3E-01 


3.7E-02 


2.4E-01 


8.6E-02 


3.5E-04 


4.4E-03 


2.1E-03 


8.5E-04 



n5b0.30vl0tl6 


1.3E+01 


1.4E+01 


2.2E+00 


2.3E+00 


9.9E+00 


7.7E+00 


6.1E+00 


2.6E+00 


1.5E-01 


8.2E-02 


3.3E-02 


1.3E-02 


n5bu.30vlOt33 


1.5E+01 


1.5E+01 


2.7E4-00 


2.9E+00 


1.1E+01 


1.0E+01 


6.9E+00 


3.1E+00 


2.2E-01 


8.6E-02 


3.4E-02 


1.3E-02 


n5b0.30vl0t67 


1.8E+01 


1.6E+01 


2.9E4-00 


3.2E+00 


1.2E+01 


1.1E+01 


7.3E+00 


3.3E+00 


2.4E-01 


9.0E-02 


3.6E-02 


1.4E-02 


n5b0.30vl2t23 


1.6E+01 


1.7E+01 


3.7E+00 


3.8E+00 


1.3E+01 


1.6E+01 


8.0E+00 


4.6E4-00 


4.3E-01 


1.3E-01 


4.7E-02 


1.7E-02 


n5b0.30vl2t48 


2.2E+01 


2.1E+01 


4.6E4-00 


4.6E+00 


1.6E+01 


1.9E+01 


9.8E+00 


5.8E4-00 


5.3E-01 


1.7E-01 


6.6E-02 


2.5E-02 


n5b0.30v!2t69 


2.4E+01 


2.1E+01 


4.4E+00 


4.4E+00 


1.6E+01 


1.8E+01 


9.9E+00 


5.6E+00 


4.8E-01 


1.8E-01 


6.9E-02 


2.7E-02 



n5b().45vl0t51 


1.4E+01 


1.4E+01 


2.5E+00 


2.7E+00 


9.9E4-00 


9.4E+00 


6.0E+00 


2.7E4-00 


2.0E-01 


7.1E-02 


2.8E- 


•02 


LIE- 


■02 


n5b0.45vl0t88 


1.5E+01 


1.4E+01 


2.6E4-00 


2.9E+00 


1.0E+01 


9.9E+00 


6.2E+00 


2.8E4-00 


2.2E-01 


7.2E-02 


2.9E- 


•02 


LIE- 


■02 


n5b0.45vl0tll5 


1.6E+01 


1.5E+01 


2.6E4-00 


2.9E+00 


1.0E+01 


1.0E+01 


6.3E+00 


2.8E4-00 


2.2E-01 


7.4E-02 


3.0E- 


■02 


LIE- 


■02 


n5b0.45vl2t58 


1.5E+01 


1.5E+01 


3.4E+00 


3.6E+00 


1.2E+01 


1.4E+01 


7.0E+00 


3.9E4-00 


3.7E-01 


1.0E-01 


3.8E- 


•02 


1.4E- 


■02 


n5b0.45vl2t76 


1.8E4-01 


1.7E+01 


3.7E+00 


3.8E+00 


1.3E+01 


1.5E+01 


7.7E+00 


4.3E4-00 


4.0E-01 


1.2E-01 


4.6E- 


■02 


1.7E- 


■02 


n5b0.45vl5t24 


1.4E+01 


1.5E+01 


3.6E4-00 


3.4E+00 


1.1E+01 


1.7E+01 


7.1E+00 


4.7E4-00 


5.0E-01 


1.5E-01 


5.0E- 


■02 


1.8E- 


■02 


n5b0.45vl5t40 


2.5E+01 


2.6E+01 


7.3E+00 


6.8E+0O 


2.1E+01 


3.6E4-01 


1.4E+01 


9.5E4-00 


1.2E+00 


3.9E-01 


1.3E- 


■01 


4.6E- 


■02 


n5b0.45vl5t55 


3.3E+01 


3.2E+01 


8.3E+00 


7.8E+00 


2.5E4-01 


4.0E4-01 


1.6E+01 


1.1E+01 


1.4E+00 


5.0E-01 


1.7E- 


•01 


6.5E 


■02 


n5b0.45vl8tl3 


1.6E+01 


1.8E+01 


4.7E4-00 


4.3E+00 


1.5E+01 


2.4E4-01 


9.6E+00 


7.4E4-00 


7.7E-01 


4.4E-01 


1.5E- 


■01 


5.4E- 


■02 


n5b0.45vl8t24 


3.5E+01 


4.2E+01 


1.4E+01 


1.2E+01 


3.4E4-01 


7.6E4-01 


2.2E+01 


1.9E+01 


2.9E4-00 


2.1E+00 


6.4E- 


■01 


2.0E- 


■01 


n5b().45v!8t35 


4.8E+01 


5.4E+01 


1.7E+01 


1.5E+01 


4.3E+01 


8.9E+01 


2.9E+01 


2.3E4-01 


4.1E+00 


2.4E+00 


6.8E- 


■01 


2. IE. 


■01 



Table B.l. (continued) Predicted integrated intensities of the 0-H2O lines to be targeted by the HIFI and PACS receivers of the Herschel telescope. The units are K km s _1 . The 
models considered are those which best fit the H2 pure rotational excitation diagram at the position of the SiO knot. A model referred to as 'n4b0.45vl5t335' has the following 
set of input parameters: a pre-shock density of 10 4 cm 3 , a magnetic-field parameter of 0.45, a shock velocity of 15 km s _1 , and an age of 335 years. 



transition name 


1 ,r. 1„, 

r 10 - J01 


3 12 " 303 


7~, 
312 " ^21 


7~, 9., 
Z21 " Z12 


7,~ 

zi2 - ioi 


303 " Z12 


J30 " J21 


4-14 - 3o3 


Z2i - r 10 


7o7 - 6l6 


G . „ 7- 
018 " '07 


yoi - 


018 


frequency (GHz) 


JJO.VJO 


1 HQ7 


|1« 1 77 
1 1 J J. 1Z/ 


lOOl.UUo 


1 AAQ QO^ 

looy.yuj 


1 71 A 77H 

1 / 10. / /u 


9 1 OA 1 Afx 


7AAH A1A 
Z04U.4 /4 


7771 Q77 
Z / / J.y 1 1 


A\ AA 8^7 
4100.0JZ 


AI'XA 7QA 
4 /74.ZV0 


^77A 


JZU 


n5b0.60vl0t52 


1.3E+01 


1.2E+01 


2.1E+00 


2.3E+00 


8.7E+00 


7.8E+00 


5.OE+00 


2.2E+00 


1.6E-01 


6.0E-02 


2.4E-02 


9.2E- 


■03 


n5b0.60vl0tl05 


1.4E+01 


1.3E+01 


2.3E+00 


2.6E+00 


9.2E+00 


8.9E+00 


5.4E+00 


2.4E+00 


1.9E-01 


6.1E-02 


2.5E-02 


9.3E 


■03 


n5b0.60vl0tl55 


1.4E+01 


1.3E+01 


2.4E+00 


2.7E+00 


9.3E+00 


9.1E+00 


5.5E+00 


2.4E+00 


1.9E-01 


6.2E-02 


2.5E-02 


9.4E- 


■03 


n5b0.60vl2t77 


1.3E+01 


1.3E+01 


2.6E+00 


2.8E+00 


9.4E+00 


1.1E+01 


5.6E+00 


2.8E+00 


2.5E-01 


7.1E-02 


2.7E-02 


1.0E- 


■02 


n5b0.60vl2t85 


1.3E+01 


1.3E+01 


2.6E+00 


2.8E+00 


9.4E+00 


1.1E+01 


5.6E+00 


2.8E+00 


2.6E-01 


7.1E-02 


2.8E-02 


1.0E- 


■02 


n5b0.60vl2tl20 


1.7E+01 


1.5E+01 


3.2E+00 


3.4E+00 


1.1E+01 


1.3E+01 


6.7E+00 


3.6E+00 


3.4E-01 


9.7E-02 


3.7E-02 


1.4E- 


■02 


n5b0.6()vl5t53 


1.5E+01 


1.5E+01 


3.9E+00 


3.8E+00 


1.1E+01 


1.9E+01 


7.4E+00 


4.8E+00 


5.5E-01 


1.4E-01 


4.5E-02 


1.6E- 


■02 


n5b0.60vl5t72 


2.5E+01 


2.5E+01 


6.7E+00 


6.4E+00 


1.9E+01 


3.2E+01 


1.3E+01 


8.4E+00 


1.1E+00 


3.2E-01 


1.1E-01 


3.9E- 


■02 


n5b().60v20t26 


UE+01 


1.3E+01 


4.0E+00 


3.4E+00 


1.0E+01 


2.3E+01 


6.6E+00 


5.9E+00 


7.5E-01 


5.9E-01 


2.0E-01 


6.8E- 


■02 


n5b0.60v20t34 


3.1E+01 


3.5E+01 


1.2E+01 


1.0E+01 


2.7E+01 


6.8E+01 


1.8E+01 


1.5E+01 


3.3E+00 


3.0E+00 


1.1E+00 


3.1E- 


•01 


n5b0.60v20t62 


6.2E+01 


7.0E+01 


2.3E+01 


2.0E+01 


5.5E+01 


1.2E+02 


3.7E+01 


2.9E+01 


6.9E+00 


4.8E+00 


1.5E+00 


4.0E- 


■01 



n5b0.75vl0tll5 


1.3E+01 


1.2E+01 


2.1E+00 


2.4E+00 


8.4E+00 


7.9E+00 


4.8E+00 


2.1E+00 


1.6E-01 


5.4E-02 


2.2E-02 


8.2E-03 


n5b0.75vl0t220 


1.3E+01 


1.2E+01 


2.2E+00 


2.5E+00 


8.5E+00 


8.4E+00 


5.0E+00 


2.2E+00 


1.8E-01 


5.4E-02 


2.2E-02 


8.2E-03 


n5b0.75vl0t250 


1.3E+01 


1.2E+01 


2.2E+00 


2.5E+00 


8.5E+00 


8.4E+00 


5.0E+00 


2.2E+00 


1.8E-01 


5.4E-02 


2.2E-02 


8.2E-03 


n5b0.75vl2tl00 


1.3E+01 


1.2E+01 


2.4E+00 


2.7E+00 


8.9E+00 


9.8E+00 


5.2E+00 


2.6E+00 


2.3E-01 


6.6E-02 


2.6E-02 


9.6E-03 


n5b0.75vl2tl40 


1.3E+01 


1.2E+01 


2.5E+00 


2.7E+00 


9.0E+00 


1.0E+01 


5.3E+00 


2.6E+00 


2.4E-01 


6.6E-02 


2.6E-02 


9.6E-03 


n5b0.75vl2tl65 


1.4E+01 


1.3E+01 


2.7E+00 


3.0E+00 


9.6E+00 


UE+01 


5.7E+00 


2.9E+00 


2.7E-01 


7.4E-02 


2.9E-02 


1.1E-02 


n5b0.75vl5t86 


1.4E+01 


1.4E+01 


3.5E+00 


3.5E+00 


1.0E+01 


1.6E+01 


6.5E+00 


4.0E+00 


4.8E-01 


1.0E-01 


3.5E-02 


1.2E-02 


n5b().75vl5tl04 


2.1E+01 


2.1E+01 


5.7E+00 


5.5E+00 


1.6E+01 


2.7E+01 


1.1E+01 


6.8E+00 


9.0E-01 


2.3E-01 


7.7E-02 


2.8E-02 



n5bl.00vl0t385 


1.3E+01 


1.1E+01 


1.9E+00 


2.3E+00 


7.6E+00 


7.5E+00 


4.4E+00 


1.8E+00 


1.5E-01 


4.5E-02 


1.8E-02 


6.7E-03 


n5bl.00vl0t410 


1.3E+01 


UE+01 


1.9E+00 


2.3E+00 


7.7E+00 


7.5E+00 


4.4E+00 


1.8E+00 


1.5E-01 


4.4E-02 


1.8E-02 


6.6E-03 


n5bl.00vl2t205 


1.3E+01 


1.2E+01 


2.3E+00 


2.6E+00 


8.5E+00 


9.3E+00 


4.9E+00 


2.4E+00 


2.1E-01 


6.0E-02 


2.4E-02 


9.0E-03 


n5bl.00vl2t235 


1.3E+01 


1.2E+01 


2.3E+00 


2.6E+00 


8.4E+00 


9.3E+00 


4.9E+00 


2.4E+00 


2.1E-01 


6.0E-02 


2.4E-02 


8.9E-03 


n5bl.00vl2t255 


1.3E+01 


1.2E+01 


2.3E+00 


2.6E+00 


8.4E+00 


9.3E+00 


4.9E+00 


2.4E+00 


2.1E-01 


6.0E-02 


2.4E-02 


8.9E-03 


n5bl.00vl5tl60 


1.3E+01 


1.2E+01 


3.1E+00 


3.2E+00 


9.2E+00 


1.4E+01 


5.8E+00 


3.4E+00 


4.0E-01 


8.0E-02 


2.8E-02 


1.0E-02 


n5bl.00vl5tl75 


1.6E+01 


1.5E+01 


3.9E+00 


4.0E+00 


UE+01 


1.8E+01 


7.3E+00 


4.4E+00 


5.4E-01 


1.2E-01 


4.0E-02 


1.4E-02 



n5bl.25vl5t205 1.2E+01 UE+01 2.6E+00 2.8E+00 8.1E+00 1.2E+01 4.9E+00 2.7E+00 3.1E-01 6.2E-02 2.2E-02 8.1E-03 
n5bl.25vl5t265 1.4E+01 1.3E+01 3.3E+00 3.5E+00 9.7E+00 1.5E+01 6.0E+00 3.5E+00 4.2E-01 8.7E-02 3.1E-02 1.1E-02 



n5bl.50vl0t525 UE+01 9.8E+00 1.6E+00 2.0E+00 6.4E+00 6.0E+00 3.5E+00 1.4E+00 UE-01 3.1E-02 1.2E-02 4.4E-03 

n5bl.50vl0t775 UE+01 9.7E+00 1.6E+00 2.0E+00 6.4E+00 6.1E+00 3.5E+00 1.4E+00 1.2E-01 3.0E-02 1.2E-02 4.4E-03 

n5bl.50vl5t335 UE+01 UE+01 2.4E+00 2.7E+00 7.7E+00 UE+01 4.6E+00 2.5E+00 2.7E-01 5.6E-02 2.1E-02 7.5E-03 

n5bl.50vl5t370 1.2E+01 UE+01 2.6E+00 2.8E+00 8.0E+00 UE+01 4.8E+00 2.6E+00 3.0E-01 6.0E-02 2.2E-02 8.0E-03 



Table B.l. (continued) Predicted integrated intensities of the 0-H2O lines to be targeted by the HIFI and PACS receivers of the Herschel telescope. The units are K km s _1 . The 
models considered are those which best fit the H2 pure rotational excitation diagram at the position of the SiO knot. A model referred to as 'n4b0.45vl5t335' has the following 
set of input parameters: a pre-shock density of 10 4 cm 3 , a magnetic-field parameter of 0.45, a shock velocity of 15 km s _1 , and an age of 335 years. 



transition name 


r 10 - J01 


J12 " J03 


3l2 - Z21 


9~, 9., 
Z21 " Z12 


9,~ 

Z12 - ioi 


J03 " Z12 


930 " J21 


4l4 - 3o3 


^21 " A 10 


7()7 - 6l6 


0I8 " '07 


yoi - sis 


frequency (GHz) 


JJO.VJO 


iuy/.joj 


1 1 1 97 
1 1 J J. 1Z / 


lOOl.UUo 


1 aaq an^ 
looy.yuj 


1 71 £ 77H 


9 1 OA 1 AA 


9AAH A1A 
Z04U.4 /4 


Li / j.y / 1 


A\ AA 8^9 




C97A ^on 
JZ /O. JZU 


n5bl.75vl0t505 


1.0E+01 


9.1E+00 


1.4E+00 


1.8E+00 


5.8E+00 


5.3E+00 


3.1E+00 


1 .2E+00 


9.8E-02 


2.5E-02 


9.9E-03 


3.5E-03 


n5bl.75vl0t915 


1.0E+01 


9.1E+00 


1.4E+00 


1.9E+00 


5.8E+00 


5.5E+00 


3.2E+00 


1.2E+00 


1.0E-01 


2.5E-02 


9.7E-03 


3.4E-03 


n5bl.75vl2t315 


1.1E+01 


1.0E+01 


2.3E+00 


2.5E+00 


7.5E+00 


9.7E+00 


4.3E+00 


2.3E+00 


2.4E-01 


5.3E-02 


2.0E-02 


7.4E-03 


n5bl.75vl2t410 


1.1E+01 


1.0E+01 


2.3E+00 


2.5E+00 


7.5E+00 


9.9E+00 


4.3E+00 


2.3E+00 


2.5E-01 


5.3E-02 


2.0E-02 


7.3E-03 


n5bl.75vl2t465 


1.1E+01 


1.0E+01 


2.3E+00 


2.6E+00 


7.5E+00 


9.9E+00 


4.4E+00 


2.3E+00 


2.5E-01 


5.3E-02 


2.0E-02 


7.3E-03 


n5bl.75vl8t305 


1.3E+01 


1.2E+01 


3.4E+00 


3.4E+00 


8.8E+00 


1.6E+01 


6.0E+00 


3.5E+00 


4.8E-01 


8.6E-02 


2.7E-02 


9.3E-03 


n5bl.75vl8t325 


1.3E+01 


1.2E+01 


3.4E+00 


3.4E+00 


8.9E+00 


1.6E+01 


6.1E+00 


3.5E+00 


4.9E-01 


8.7E-02 


2.8E-02 


9.4E-03 


n5bl.75v!8t340 


1.4E+01 


1.3E+01 


3.8E+00 


3.7E+00 


9.6E+00 


1.8E+01 


6.6E+00 


3.9E+00 


5.7E-01 


1.0E-01 


3.2E-02 


1.1E-02 



n5b2.00vl8t435 1.2E+01 1.1E+01 3.0E+00 3.1E+00 8.1E+00 1.4E+01 5.4E+00 3.0E+00 4.0E-01 7.0E-02 2.3E-02 8.0E-03 
n5b2.00v!8t440 1.2E+01 UE+01 3.1E+00 3.2E+00 8.3E+00 1.4E+01 5.5E+00 3.1E+00 4.2E-01 7.3E-02 2.4E-02 8.4E-03 



A. Gusdorf et al.: Revisiting the shocks in BHR71: new observational constraints and H 2 predictions for Herschel 

Table B.2. Predicted integrated intensities of the P-H2O lines to be targeted by the HIFI and PACS receivers of the Herschel 
telescope. The units are K km s _1 . The models considered are those which best fit the H2 pure rotational excitation diagram at the 
position of the SiO knot. A model referred to as 'n4b0.45vl5t335' has the following set of input parameters: a pre-shock density of 
10 4 crrT 3 , a magnetic-field parameter of 0.45, a shock velocity of 15 km s _1 , and an age of 335 years. 



transition name 


In - Ooo 


2o2 - 111 


2n - 2o2 


3 13 " 2()2 


4o4 - 3l3 


322 - 2n 


frequency (GHz) 


1113.343 


987.927 


752.033 


2164.132 


2391.573 


3331.458 


n4b().45vl5t335 


3.3E+00 


4.2E+00 


3.4E+00 


5.4E-01 


1.9E-02 


1.4E-01 


n4b0.45vl5t555 


4.1E+00 


5.1E+00 


3.9E+00 


7.1E-01 


2.6E-02 


2.0E-01 


n4b0.45vl5t645 


4.5E+00 


5.5E+00 


4.2E+00 


7.6E-01 


2.8E-02 


2.1E-01 


n4b0.45vl8t350 


4.9E+00 


6.6E+00 


5.1E+00 


1.2E+00 


6.9E-02 


3.7E-01 


n4b0.45vl8t385 


5.1E+00 


6.7E+00 


5.2E+00 


1.2E+00 


7.1E-02 


3.8E-01 


n4b0.45vl8t425 


5.2E+00 


6.9E+00 


5.3E+00 


1.3E+00 


7.5E-02 


4.0E-01 


n4b0.45v20t210 


5.2E+00 


7.2E+00 


5.7E+00 


1.5E+00 


1.2E-01 


5.1E-01 


n4b0.45v20t315 


5.8E+00 


7.8E+00 


6.1E+00 


1.7E+00 


1.3E-01 


5.8E-01 


n4b0.45v20t600 


6.5E+00 


8.2E+00 


6.1E+00 


1.6E+00 


1.2E-01 


5.7E-01 



n4b() 


,60vl5t640 


3.6E+00 


4.4E+00 


3.4E+00 


5.6E-01 


1.8E- 


■02 


1.5E- 


■01 


n4b0 


60vl5t655 


3.6E+00 


4.5E+00 


3.4E+00 


5.7E-01 


1.9E- 


■02 


1.5E 


-01 


n4b0 


.60vl5t885 


4.4E+00 


5.1E+00 


3.8E+00 


6.7E-01 


2.2E- 


■02 


1.9E- 


-01 


n4b0 


.60v 181370 


4.3E+00 


5.7E+00 


4.5E+00 


9.8E-01 


5.2E- 


■02 


2.9E- 


01 


n4b0 


.60v 181515 


5.0E+00 


6.4E+00 


4.9E+00 


1.1E+00 


6.0E- 


■02 


3.4E- 


01 


n4b0 


60v 181730 


6.0E+00 


7.2E+00 


5.3E+00 


1.3E+00 


7.3E- 


■02 


4.1E- 


01 


n4b0 


.60v20t260 


4.8E+00 


6.5E+00 


5.2E+00 


1.3E+00 


1.0E- 


01 


4.5E- 


01 


n4b0 


.60v20t330 


5.3E+00 


6.9E+00 


5.6E+00 


1.4E+00 


LIE- 


01 


4.9E- 


01 


n4b0 


.60v20t560 


6.4E+00 


7.9E+00 


6.0E+00 


1.5E+00 


LIE- 


01 


5.4E- 


01 


n4b0 


60v22tl95 


6.0E+00 


8.4E+00 


7.1E+00 


2.0E+00 


2.3E- 


01 


8.3E- 


01 


n4b0 


.60v22t275 


6.8E+00 


9.0E+00 


7.6E+00 


2.2E+00 


2.5E- 


01 


9.0E- 


01 


n4b0 


,60v22t365 


7.0E+00 


9.1E+00 


7.5E+00 


2.1E+00 


2.3E- 


01 


8.7E- 


01 



n4b().75vl5t710 


2.9E+00 


3 


.6E+00 


2.8E+00 


4.1E-01 


1.2E- 


■02 


LIE- 


-01 


n4b0.75vl5t830 


3.3E+00 


4 


.0E+00 


3.1E+00 


4.7E-01 


1.4E- 


■02 


1.2E- 


-01 


n4b0.75vl5tll20 


3.9E+00 


4 


.4E+00 


3.2E+00 


5.2E-01 


1.5E- 


■02 


1.4E- 


01 


n4b0.75vl8t375 


3.2E+00 


4 


.3E+00 


3.5E+00 


6.8E-01 


3.1E- 


■02 


1.8E- 


01 


n4b0.75vl8t510 


4.0E+00 


5 


.2E+00 


4.1E+00 


8.5E-01 


4.1E- 


■02 


2.4E- 


01 


n4b0.75vl8t695 


5.0E+00 


6 


.2E+00 


4.7E+00 


1.0E+00 


5.0E- 


■02 


3.1E- 


01 


n4b0.75v20t365 


4.7E+00 


6 


.1E+00 


5.0E+00 


1.2E+00 


8.9E- 


■02 


4.0E- 


01 


n4b0.75v20t430 


5.0E+00 


6 


.5E+00 


5.2E+00 


1.3E+00 


9.5E- 


■02 


4.3E- 


01 


n4b0.75v20t530 


5.6E+00 


7 


.0E+00 


5.5E+00 


1.4E+00 


1.0E- 


01 


4.7E- 


01 


n4b0.75v22t225 


5.0E+00 


7 


.0E+00 


6.0E+00 


1.6E+00 


1.8E- 


01 


6.6E- 


01 


n4b0.75v22t285 


5.8E+00 


7 


.8E+00 


6.7E+00 


1.9E+00 


2.1E- 


01 


7.6E- 


01 


n4b().75v22t450 


6.9E+00 


8 


.6E+00 


7.1E+00 


2.0E+00 


2.1E- 


01 


7.8E- 


01 



n4bl.00v2()t490 


3.5E+00 


4.7E+00 


4.0E+00 


8.5E-01 


5.2E-02 


2.6E-01 


n4bl.00v20t605 


4.6E+00 


5.8E+00 


4.7E+00 


1.1E+00 


6.9E-02 


3.4E-01 


n4bl.00v20t740 


5.3E+00 


6.4E+00 


5.0E+00 


1.2E+00 


7.8E-02 


3.9E-01 


n4bl.00v22t370 


4.2E+00 


5.8E+00 


5.0E+00 


1.3E+00 


1.2E-01 


4.8E-01 


n4bl.00v22t435 


5.1E+00 


6.6E+00 


5.8E+00 


1.5E+00 


1.6E-01 


5.8E-01 


n4bl.00v22t465 


5.3E+00 


6.8E+00 


5.9E+00 


1.5E+00 


1.6E-01 


6.0E-01 


n4bl.00v25t290 


5.2E+00 


7.5E+00 


6.4E+00 


1.9E+00 


2.3E-01 


8.3E-01 


n4bl.00v25t320 


6.5E+00 


8.8E+00 


8.2E+00 


2.5E+00 


3.6E-01 


1.2E+00 


n4bl.00v25t345 


7.5E+00 


9.8E+00 


9.5E+00 


2.8E+00 


4.6E-01 


1.5E+00 



n4bl.25v22t585 


3.9E+00 


5.1E+00 


4.5E+00 


1.0E+00 


8.9E-02 


3.7E-01 


n4bl.25v22t625 


4.4E+00 


5.7E+00 


4.9E+00 


1.2E+00 


1.0E-01 


4.2E-01 


n4bl.25v22t670 


4.7E+00 


5.9E+00 


5.1E+00 


1.2E+00 


1.1E-01 


4.5E-01 


n4bl.25v25t425 


5.5E+00 


7.4E+00 


7.2E+00 


2.0E+00 


3.1E-01 


1.0E+00 


n4bl.25v25t445 


6.0E+00 


7.8E+00 


7.7E+00 


2.2E+00 


3.4E-01 


1.1E+00 


n4bl.25v25t530 


7.8E+00 


9.5E+00 


9.6E+00 


2.7E+00 


4.6E-01 


1.5E+00 


n4bl.25v28t325 


2.3E-01 


1.6E-01 


5.3E-02 


2.0E-02 


2.4E-04 


8.1E-03 


n4bl.25v28t465 


2.2E+00 


1.7E+00 


7.5E-01 


1.9E-01 


3.2E-03 


7.5E-02 
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A. Gusdorf et al.: Revisiting the shocks in BHR71: new observational constraints and H 2 predictions for Herschel 

Table B.2. (continued) Predicted integrated intensities of the P-H2O lines to be targeted by the HIFI and PACS receivers of the 
Herschel telescope. The units are K km s _1 . The models considered are those which best fit the H2 pure rotational excitation 
diagram at the position of the SiO knot. A model referred to as 'n4b0.45vl5t335' has the following set of input parameters: a 
pre-shock density of 10 4 ctrT 3 , a magnetic-field parameter of 0.45, a shock velocity of 15 km s _1 , and an age of 335 years. 



transition name 


In - Ooo 


2o2 - 111 


2n - 2o2 


3 13 " 2()2 


4()4 - 3l3 


322 - 2n 


frequency (GHz) 


1113.343 


987.927 


752.033 


2164.132 


2391.573 


3331.458 


n4bl.50v22t830 


3.3E+00 


4.4E+00 


3.9E+00 


8.2E-01 


5.8E-02 


2.6E-01 


n4bl.50v22t930 


4.3E+00 


5.4E+00 


4.6E+00 


1.0E+00 


7.7E-02 


3.5E-01 


n4bl.50v22tl005 


5.0E+00 


5.9E+00 


4.9E+00 


UE+00 


8.1E-02 


3.8E-01 


n4bl.50v25t605 


5.3E+00 


6.9E+00 


6.9E+00 


1.8E+00 


2.7E-01 


9.1E-01 


n4bl.50v25t655 


6.0E+00 


7.5E+00 


7.6E+00 


2.0E+00 


3.1E-01 


1.0E+00 


n4bl.50v25t755 


7.3E+00 


8.6E+00 


8.3E+00 


2.2E+00 


3.4E-01 


1.1E+00 


n4bl.50v28t500 


1.1E+00 


9.7E-01 


5.4E-01 


8.9E-02 


2.0E-03 


3.2E-02 


n4bl.50v28t520 


1.2E+00 


9.8E-01 


5.2E-01 


9.3E-02 


1.9E-03 


3.4E-02 


n4bl.50v28t625 


2.2E+00 


1.7E+00 


7.4E-01 


1.8E-01 


2.9E-03 


7.2E-02 


n4bl.50v30t405 


7.5E-01 


5.9E-01 


2.9E-01 


5.7E-02 


1.1E-03 


2.1E-02 


n4bl.50v30t435 


9.1E-01 


6.8E-01 


2.9E-01 


7.2E-02 


1.1E-03 


2.8E-02 


n4bl.50v30t825 


6.4E+00 


5.1E+00 


2.2E+00 


5.8E-01 


1.0E-02 


2.4E-01 



n4bl.75v22tlll5 


3.0E+00 


3.9E+00 


3.4E+00 


6.9E-01 


4.5E- 


■02 


2.2E-01 


n4bl.75v22tll80 


3.5E+00 


4.4E+00 


3.8E+00 


7.7E-01 


4.8E- 


■02 


2.4E-01 


n4bl.75v22tl260 


4.1E+00 


5.0E+00 


4.1E+00 


8.7E-01 


5.6E- 


■02 


2.8E-01 


n4bl.75v25t855 


5.2E+00 


6.5E+00 


6.5E+00 


1.6E+00 


2.3E- 


01 


7.8E-01 


n4bl.75v25t885 


5.4E+00 


6.7E+00 


6.7E+00 


1.7E+00 


2.4E- 


01 


8.1E-01 


n4bl.75v25t920 


5.8E+00 


7.1E+00 


6.9E+00 


1.8E+00 


2.4E- 


01 


8.3E-01 


n4bl.75v28t660 


9.7E-01 


9.1E-01 


5.9E-01 


7.7E-02 


2.1E- 


■03 


2.5E-02 


n4bl.75v28t7()5 


1.3E+00 


1.1E+00 


5.9E-01 


1.0E-01 


2.1E- 


■03 


3.7E-02 


n4bl.75v28t890 


9.6E+00 


1.1E+01 


1.0E+01 


2.8E+00 


4.2E- 


01 


1.4E+00 


n4bl.75v30t570 


8.1E-01 


6.3E-01 


3.1E-01 


6.0E-02 


1.0E- 


■03 


2.2E-02 


n4bl.75v30t600 


1.1E+00 


8.4E-01 


3.7E-01 


8.7E-02 


1.3E- 


■03 


3.3E-02 


n4bl.75v3()t670 


2.0E+00 


1.5E+00 


6.4E-01 


1.7E-01 


2.6E- 


■03 


6.7E-02 



n4b2.()0v28t835 


1.1E+00 


1.2E+00 


9.3E-01 


1.1E-01 


3.6E-03 


3.2E-02 


n4b2.00v28t880 


1.5E+00 


1.6E+00 


1.2E+00 


1.6E-01 


4.6E-03 


4.8E-02 


n4b2.00v28t920 


3.1E+00 


3.7E+00 


2.8E+00 


4.9E-01 


1.7E-02 


1.4E-01 


n4b2.00v30t730 


7.2E-01 


5.8E-01 


3.2E-01 


5.0E-02 


9.7E-04 


1.7E-02 


n4b2.00v30t755 


9.5E-01 


7.7E-01 


4.0E-01 


7.1E-02 


1.3E-03 


2.6E-02 


n4b2.()0v3()t775 


1.1E+00 


8.6E-01 


4.0E-01 


8.6E-02 


1.4E-03 


3.2E-02 



n5b0.30vl0tl6 


1.1E+01 


1.5E+01 


1.1E+01 


2.8E+00 


1.9E-01 


9.9E-01 


n5b0.30vlOt33 


1.2E+01 


1.7E+01 


1.3E+01 


3.4E+00 


2.4E-01 


1.2E+00 


n5b0.30vl0t67 


1.3E+01 


1.8E+01 


1.4E+01 


3.7E+00 


2.7E-01 


1.3E+00 


n5b0.30vl2t23 


1.3E+01 


1.9E+01 


1.4E+01 


4.7E+00 


4.6E-01 


1.8E+00 


n5b0.30vl2t48 


1.6E+01 


2.3E+01 


1.7E+01 


5.8E+00 


5.8E-01 


2.3E+00 


n5b().30vl2t69 


1.7E+01 


2.3E+01 


1.7E+01 


5.6E+00 


5.3E-01 


2.2E+00 



n5b().45vl0t51 


1.1E+01 


1.6E+01 


1.2E+01 


3.0E+00 


2.0E-01 


9.9E-01 


n5b0.45vl0t88 


1.2E+01 


1.6E+01 


1.2E+01 


3.2E+00 


2.2E-01 


1.0E+00 


n5b0.45vl0tll5 


1.2E+01 


1.6E+01 


1.3E+01 


3.2E+00 


2.2E-01 


1.1E+00 


n5b0.45vl2t58 


1.2E+01 


1.7E+01 


1.3E+01 


4.1E+00 


3.8E-01 


1.5E+00 


n5b0.45vl2t76 


1.3E+01 


1.9E+01 


1.4E+01 


4.5E+00 


4.2E-01 


1.7E+00 


n5b0.45vl5t24 


1.1E+01 


1.6E+01 


1.2E+01 


4.5E+00 


6.0E-01 


2.0E+00 


n5b0.45vl5t40 


2.0E+01 


2.8E+01 


2.2E+01 


9.0E+00 


1.6E+00 


4.7E+00 


n5b0.45vl5t55 


2.4E+01 


3.3E+01 


2.5E+01 


1.0E+01 


1.8E+00 


5.4E+00 


n5b0.45vl8tl3 


1.3E+01 


2.0E+01 


1.5E+01 


6.3E+00 


1.0E+00 


3.2E+00 


n5b0.45vl8t24 


2.8E+01 


4.3E+01 


3.4E+01 


1.7E+01 


4.1E+00 


1.1E+01 


n5b0.45vl8t35 


3.8E+01 


5.4E+01 


4.3E+01 


2.1E+01 


5.3E+00 


1.3E+01 
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Table B.2. (continued) Predicted integrated intensities of the P-H2O lines to be targeted by the HIFI and PACS receivers of the 
Herschel telescope. The units are K km s _1 . The models considered are those which best fit the H2 pure rotational excitation 
diagram at the position of the SiO knot. A model referred to as 'n4b0.45vl5t335' has the following set of input parameters: a 
pre-shock density of 10 4 ctrT 3 , a magnetic-field parameter of 0.45, a shock velocity of 15 km s _1 , and an age of 335 years. 



transition name 


In - Ooo 


2o2 - 111 


2n - 2o2 


3 13 - 2o2 


4o4 - 3l3 


322 - 2n 


frequency (GHz) 


1113.343 


987.927 


752.033 


2164.132 


2391.573 


3331.458 


n5b0.60vl0t52 


1.0E+01 


1.4E+01 


1.1E+01 


2.5E+00 


1.5E-01 


8.1E-01 


n5b0.60vl0tl05 


1.1E+01 


1.5E+01 


1.1E+01 


2.8E+00 


1.8E-01 


8.9E-01 


n5b().60vl0tl55 


1.1E+01 


1.5E+01 


1.2E+01 


2.9E+00 


1.8E-01 


9.1E-01 


n5b0.6()vl2t77 


1.0E+01 


1.4E+01 


1.1E+01 


3.1E+00 


2.4E-01 


1.1E+00 


n5b0.60vl2t85 


1.0E+01 


1.4E+01 


1.1E+01 


3.2E+00 


2.5E-01 


1.1E+00 


n5b0.60vl2tl20 


1.2E+01 


1.7E+01 


1.3E+01 


3.9E+00 


3.4E-01 


1.4E+00 


n5b0.60vl5t53 


1.1E+01 


1.6E+01 


1.3E+01 


4.7E+00 


6.8E-01 


2.1E+00 


n5b0.60vl5f72 


1.9E+01 


2.6E+01 


2.1E+01 


8.1E+00 


1.5E+00 


4.2E+00 


n5b0.60v20t26 


8.4E+00 


1.4E+01 


1.1E+01 


5.1E+00 


1.1E+00 


3.0E+00 


n5b0.60v20t34 


2.4E+01 


3.4E+01 


2.8E+01 


1.4E+01 


4.4E+00 


1.0E+01 


n5b0.60v20t62 


5.1E+01 


6.9E+01 


5.5E+01 


2.7E+01 


8.3E+00 


1.9E+01 



n5b().75vl0tll5 


1.0E+01 


1.4E+01 


1.1E+01 


2.5E+00 


1.5E-01 


7.7E-01 


n5b0.75vl0t220 


1.0E+01 


1.4E+01 


1.1E+01 


2.6E+00 


1.6E-01 


8.1E-01 


n5b0.75vl0t250 


1.0E+01 


1.4E+01 


1.1E+01 


2.6E+00 


1.6E-01 


8.1E-01 


n5b0.75vl2tl00 


9.7E+00 


1.4E+01 


1.1E+01 


2.9E+00 


2.1E-01 


9.7E-01 


n5b0.75vl2tl40 


9.8E+00 


1.4E+01 


1.1E+01 


3.0E+00 


2.2E-01 


1.0E+00 


n5b0.75vl2tl65 


1.0E+01 


1.5E+01 


1.1E+01 


3.2E+00 


2.5E-01 


1.1E+00 


n5b0.75vl5t86 


1.0E+01 


1.5E+01 


1.2E+01 


4.1E+00 


5.6E-01 


1.8E+00 


n5b().75vl5tl04 


1.6E+01 


2.2E+01 


1.8E+01 


6.7E+00 


1.1E+00 


3.3E+00 



n5bl.00vl0t385 


9.4E+00 


1.3E+01 


1.0E+01 


2.3E+00 


1.3E-01 


6.8E-01 


n5bl.00vl0t410 


9.4E+00 


1.3E+01 


1.0E+01 


2.3E+00 


1.3E-01 


6.8E-01 


n5bl.00vl2t205 


9.6E+00 


1.4E+01 


1.1E+01 


2.7E+00 


1.9E-01 


8.9E-01 


n5bl.00vl2t235 


9.5E+00 


1.3E+01 


1.0E+01 


2.7E+00 


1.9E-01 


8.9E-01 


n5bl.00vl2t255 


9.6E+00 


1.4E+01 


1.0E+01 


2.7E+00 


1.9E-01 


9.0E-01 


n5bl.00vl5tl60 


9.6E+00 


1.4E+01 


1.1E+01 


3.6E+00 


4.4E-01 


1.5E+00 


n5bl.00vl5tl75 


1.2E+01 


1.6E+01 


1.3E+01 


4.6E+00 


6.4E-01 


2.0E+00 



n5bl.25vl5t205 8.7E+00 1.2E+01 9.8E+00 3.0E+00 3.0E-01 1.1E+00 
n5bl.25vl5t265 1.0E+01 1.4E+01 1.2E+01 3.8E+00 4.5E-01 1.5E+00 



n5bl.50vl()t525 8.4E+00 1.1E+01 9.2E+00 1.8E+00 8.6E-02 5.0E-01 

n5bl.50vl0t775 8.4E+00 1.1E+01 9.2E+00 1.8E+00 8.8E-02 5.0E-01 

n5bl.50vl5t335 8.4E+00 1.2E+01 9.5E+00 2.8E+00 2.6E-01 1.0E+00 

n5bl.50vl5t370 8.7E+00 1.2E+01 9.9E+00 3.0E+00 2.8E-01 1.1E+00 



n5bl.75vl0t505 7 

n5bl.75vl0t915 7 

n5bl.75vl2t315 8 

n5bl.75vl2t410 8 

n5bl.75vl2t465 8 

n5bl.75vl8t305 9 

n5bl.75vl8t325 9 

n5bl.75vl8t340 1 



.9E+00 


1.1E+01 


.9E+00 


1.1E+01 


.3E+00 


1.2E+01 


.3E+00 


1.2E+01 


.3E+00 


1.2E+01 


.4E+00 


1.3E+01 


.5E+00 


1.3E+01 


.0E+01 


1.4E+01 



8.7E+00 


1.6E+00 


8.8E+O0 


1.6E+00 


9.3E+00 


2.6E+00 


9.4E+00 


2.6E+00 


9.4E+00 


2.6E+00 


1.1E+01 


3.7E+00 


1.1E+01 


3.7E+00 


1.2E+01 


4.1E+00 



7.0E-02 


4.3E-01 


7.3E-02 


4.4E-01 


2.1E-01 


8.9E-01 


2.2E-01 


9.1E-01 


2.2E-01 


9.2E-01 


5.7E-01 


1.7E+00 


5.8E-01 


1.8E+00 


6.8E-01 


2.0E+00 



n5b2.00vl8t435 8.8E+00 1.2E+01 1.0E+01 3.3E+00 4.5E-01 1.4E+00 
n5b2.()0vl8t440 9.0E+00 1.2E+01 1.1E+01 3.4E+00 4.7E-01 1.5E+00 
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